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PREFACE 

This  work  is  essentially  an  exposition  of  the  simplex 
algorithm  as  used  in  hand-computed  solutions  of  linear 
programming  problems.  As  originally  developed  by 
George  B.  Dantzig,  this  method  was  publicly  presented  in 
Monograph  No.  13  of  the  Cowles  Commission  for  Research 
in  Economics  in  1951.  Later  the  applicability  of  the  tech- 
nique to  industrial  problems  was  demonstrated  in  several 
original  instances  by  A.  Charnes,  W.  W.  Cooper,  and  their 
associates.  From  these  studies  an  expanding  interest  in 
the  general  usefulness  of  the  device  rapidly  developed. 
This  interest  was  augmented  by  the  publication  in  1953  by 
Charnes  et  at.  oi  An  Introduction  to  Linear  Programming. 

For  some  time  this  text  remained  the  only  exposition 
available  that  was  not  precisely  formulated  in  tight  mathe- 
matical terms  in  its  entirety.  Indeed,  the  latter  half  of 
the  work  presented  a  rigorous  mathematical  proof  of  the 
correctness  of  the  technique.  This  proof  uses  the  concepts 
and  notations  of  modern  algebra  which  are  as  yet  not 
commonly  taught  to  undergraduate  engineering  and 
science  students,  and,  since  these  concepts  have  been  de- 
veloped within  the  past  fifteen  years,  they  are  also  largely 
unknown  to  those  in  industry  and  management  whose 
college  degrees  were  taken  within  the  past  decade. 

The  present  text  was  written,  therefore,  to  explain 
Dantzig's  simplex  technique  to  persons  possessing  some 
mathematical  background — perhaps  at  the  level  of  college 
algebra.  This  group  would  comprise  all  engineers  and 
many  students  in  other  fields  of  the  arts  and  sciences. 
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The  attempt  has  been  made  to  explain  (not  prove)  for 
this  audience  how  the  technique  works  and  to  illustrate 
the  types  of  problems  to  which  it  has  been  applied.  The 
graduate  mathematician,  especially  if  his  training  has 
included  work  in  modern  algebra,  will  find  this  book  very 
elementary  and  may  be  amazed  at  the  slow  development 
of  the  concepts  as  well  as  at  their  expanded  phrasing.  He 
will  be  much  happier  with  other  presentations,  such  as  the 
one  given  in  Charnes'  book  or  perhaps  with  another  by 
Eisemann. 

The  present  volume  contains  a  number  of  sample 
problems  graded  in  difficulty  and  so  introduced  as  to 
entice  the  reader  to  try  his  hand  at  the  technique  of  prob- 
lem formulation  and  solution,  and  thus  check  his  under- 
standing of  the  ideas  involved.  This  material  has  been 
used  with  several  groups  of  first-year  graduate  students  in 
engineering  and  was  found  by  them  to  be  largely  self- 
teaching.  The  author  would  not  hesitate  to  use  it  at  the 
junior  level  with  engineering  and  science  students. 

It  is  anticipated  that  a  wider  understanding  of  this 
algorithm  will  permit  many  to  formulate  and  solve  prob- 
lems which  they  currently  believe  insoluble.  Formulation 
of  problems,  alone,  in  linear  form  may  help  to  give  many  a 
clearer  understanding  of  complex  relationships  and  a 
better  grasp  of  what  a  particular  problem  comprises. 
The  writer  is  aware  that  many  problems  can  be  formu- 
lated which  the  simplex,  even  when  machine  programmed, 
cannot  reasonably  solve  from  a  time  standpoint;  he  is 
similarly  aware  that  many  industrial  problems  are  non- 
linear. But  this  will  become  clear  also  to  the  student  who 
perseveres  and  finishes  the  problems  presented  herein. 

This  text  will  have  served  its  purpose  if  it  helps  to  bring 
this  particular  device  from  the  realm  of  mathematical 
magic  (where  some  would  place  it) ;    induce  any  to  attempt 
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more  precise  formulation  of  many-faceted  engineering, 
business,  and  management  problems;  or  even  induce  any 
to  accept  more  readily  the  aid  of  mathematicians  in  prob- 
lems beyond  the  scope  of  this  algorithm. 

The  author  wishes  to  acknowledge  with  sincere  thanks 
the  comments  and  many  helpful  criticisms  given  by  his 
students  to  whom  this  presentation  was  first  made.  As  a 
result  of  their  help,  many  unhappy  and  unclear  expressions 
were  improved.  The  manuscript  was  carefully  and  rapidly 
typed  by  Miss  Irene  Smura,  to  whom  I  also  express  my 
sincere  appreciation. 

D.  U.  Greenwald 

Syracuse,  New  York 
February,  1957 
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LINEAR  PROGRAMMING 

Introduction 

During  the  past  decade,  there  has  arisen  a  considerable 
interest  in  the  application  of  new  mathematical  methods 
to  problems  in  business  engineering  and  management, 
together  with  the  application  of  older  mathematical 
procedures  to  new  problems.  There  has  been  an  increased 
use  of  probability  in  developing  statistics  for  many  engi- 
neering applications,  with  particular  emphasis  on  electronic 
communication  and  on  statistical  quality  control.  Infor- 
mation theory  has  been  developed  as  a  procedure  for 
studying  the  transmission  of  intelligence.  Also  worthy  of 
mention  are  some  tentative  applications  of  queueing 
theory  and  game  theory.  Boolean  algebra  has  been  re- 
discovered, and  study  of  its  use  in  engineering  and  manage- 
ment is  being  urged.  And  lastly,  there  is  the  considerable 
development  of  matrix  algebra,  which  has  proved  useful 
in  electrical  engineering  and  which  gives  promise  of  being 
useful  in  dealing  with  some  other  industrial  areas.  A 
derivative  of  matrix  algebra  called  linear  programming 
appears  most  encouraging  for  these  uses,  and  it  is  to  this 
subject  that  this  book  is  devoted. 

Linear  programming  is  defined  as  the  process  of  deter- 
mining a  program  of  activity  by  finding  the  optimum 
solution  of  a  group  of  restrictive  linear  equations.  It  is 
further  required  that  the  optimum  condition  also  be 
expressed  in  equation  form.  The  conditions  of  the  prob- 
lem to  be  solved  by  this  method  must  be  such  that  they 
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can  be  reduced  to  linear  equation  form.  The  unknowns 
in  the  equations,  when  finally  determined,  will  constitute 
the  program.  The  desirability  of  the  program  will  be  a 
resultant  of  these  unknowns  and  will  be  tested  by  trying 
them  in  the  equation  which  states  the  optimum  desired. 

Less  specifically,  linear  programming  may  be  described 
as  a  mathematical  technique  for  determining  the  most 
effective,  desirable,  or  profitable  course  of  action  to  take 
when  the  situation  is  governed  by  many  known  variables 
and  conditions.  The  thought  now  occurs,  "What  can  it 
be  used  for  in  engineering  and  management?"  This 
question  has  not  as  yet  been  fully  answered.  Linear 
programming  has  been  used  to  determine  the  most  profit- 
able mixes  of  raw  materials  and  products,  notably  in  the 
petroleum  industry.*  In  simple  cases,  it  has  been  used  in 
scheduling  work  to  machines — in  machine  routing.  It 
appears  possible  to  use  it  in  inventory  control,  in  the 
location  of  warehousing,  and  in  shipping  operations.  It 
has  also  been  used,  with  some  mathematical  refinements, 
in  wage  policy  determinations.  In  all  of  these  instances, 
it  has  given  the  assurance  that  the  best  possible  answer 
was  achieved.  But  as  an  applied  tool,  it  is  still  very  close 
to  its  birth  in  the  research  situation. 

Several  computation  methods,  or  algorithms,  have  been 
developed  for  the  solution  of  linear  programming  problems. 
The  one  which  appears  most  useful  to  the  engineer  or 
manager  as  a  technique  for  exploring  the  applicability 
of  this  subject  to  the  problems  with  which  he  is  faced  is 
the   "simplex  tableau." f     This   particular  algorithm   is 

*  Reference  3. 

f  The  "simplex"  technique  was  first  presented  by  G.  B.  Dantzig  (Refer- 
ence 1).  It  was  further  modified  by  A.  Orden,  A.  J.  Hoffman,  and 
A.  Charnes.  The  latter  developed  the  mathematical  theory  in  Part  II  of 
Reference  2. 
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lOXi  +    5X2  ^  50  (d) 

6X1  +    6X2  ^  36  (e) 

4.5Xx  +  18X2  ^  81  (g) 

We  should  also  state  in  terms  of  our  unknowns  the  profit 
resulting  from  any  selected  program.  This  profit  (z) 
can  be  represented  in  dollars  as: 

9Xx  +    7X2  =  z  (/) 

This  equation  is  customarily  called  the  functional  and 
represents  the  value  we  desire  to  maximize  in  the  present 
problem. 

A  complete  and  precise  solution  of  the  example  chosen 
is  not  difficult,  since  we  have  three  simultaneous  equations 
in  only  two  unknowns.  We  need  to  make  only  a  trial  or 
two  and  check  the  program  that  maximizes  the  functional, 
i.e.,  yields  the  most  profit.  Considering  Eqs.  (a)  and  (b), 
and  solving  simultaneously,  yields  the  result  Xx  =  4, 
X2  =  2.  This  would  be  interpreted  as  a  program  of  four 
lots  of  Type  1  pins  and  two  lots  of  Type  2  pins.  But  be- 
fore adopting  this  program,  it  is  necessary  to  check  and 
see  that  limitation  (c)  is  also  satisfied.  Substituting  the 
above  lambda  values  in  (c)  gives  a  total  grinder  time  of 
54  hours,  which  is  less  than  the  limit  of  81  spare  hours 
postulated  in  the  problem.  This  substitution  also  shows 
us  that  these  three  equations  cannot  be  satisfied  simul- 
taneously— that  there  are  several  solutions.  (If  there 
were  one  solution  only,  the  right-hand  member  of  Eq.  (c) 
would  have  been  54.)  We  next  substitute  these  values  in 
our  functional  (/)  to  see  how  profitable  this  program  is. 
We  find  the  total  z  to  be  $50.00.  This  is  a  greater  profit 
than  that  obtained  from  making  Type  1  pins  (five  lots) 
only,  which  casual  inspection  of  the  problem  at  first 
suggested. 
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This  result  leads  us  to  wonder  if  some  other  values  of  X, 
which  is  to  say  some  other  program,  would  yield  a  greater 
profit.  Combining  Eqs.  (b)  and  (c)  simultaneously  gives 
Xi  as  2  lots  and  X2  as  4  lots.  Substitution  of  these  in  the 
relationship  expressed  in  Eq.  (a)  gives  a  total  right-hand 
value  of  40  hours,  which  is  less  than  the  limit  of  50  spare 
hours  postulated.  These  values  are,  then,  a  reasonable 
solution.  Checking  for  profit  with  this  program  gives  z 
as  equal  to  $36.00  which,  being  less  than  the  previous 
program,  is  less  satisfactory. 

The  last  possibility  of  simultaneous  solution  is  the 
combination  of  Eqs.  (a)  and  (c).  This  yields  program 
values  of  Xi  =  3.14  and  X2  =  3.71,  a  program  of  fractional 
lots.  Trying  these  values  in  the  relationship  expressed  in 
Eq.  (b)  gives  a  right-hand  value  of  41.1  hours.  This 
exceeds  the  limitation  of  36  hours  available  in  Lathe  B  and 
is  therefore  an  unworkable  program.  If  it  had  been  usable, 
it  would  have  yielded  a  profit  of  $54.23 — but  this  is  wishful 
thinking. 

A  clearer  picture  of  what  is  involved  in  these  limitations 
can  be  gained  by  a  plot  of  the  three  equations.  Since 
there  are  only  two  unknowns,  we  let  Xx  serve  as  the  ab- 
scissa and  X2  as  the  ordinate.  Fig.  1  shows  the  three 
straight  lines  which  are  Eqs.  (a),  (6),  and  (c).  But  it 
shows  more  than  this.  It  is  also  a  representation  of  the 
inequalities  (d),  (e),  and  (g),  if  we  consider  any  point  in 
an  area  bounded  by  an  equation  and  the  axes  as  satisfying 
the  inequation  conditions.  With  this  concept  in  mind, 
we  readily  see  that  any  pair  of  values  which  will  satisfy 
our  limiting  conditions  must  lie  within  the  area  OABCD. 
Another  way  in  which  this  fact  might  be  stated  would  be 
to  say  that  any  point  in  this  (two-dimensional)  space, 
bounded  by  the  limiting  equations  and  the  axes,  is  a  work- 
able program.    In  so  stating,  we  tie  in  our  new  concept  of  a 
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program  with  the  well-known  point  concept  in  analytic 
geometry.  And  thus  another  relationship  comes  to  view. 
The  coordinates  of  the  geometric  space  are  the  program 
values  which  satisfy  the  constraining  inequations. 


0  5  10  15  20  *, 

Fig.  1.    A  plot  of  the  limiting  conditions  for  Problem  1. 


On  the  basis  of  this  thinking,  then,  any  point  in  the 
space  (area  in  this  simple  problem)  OABCD  will  be  a 
workable  program.  Our  problem  demands  that  we  choose 
the  workable  program  which  yields  the  maximum  profit. 
This  can  be  paraphrased  neatly  to  say  that  the  solution 
must  maximize  the  functional.  To  say  it  still  another  way, 
we  want  to  choose  the  point  P  which  has  coordinates  Xx 
and  X2  such  that  z  of  Eq.  (/),  which  is  a  function  of  Xx,  X2, 
becomes  as  large  as  possible.  Our  previous  solutions  by 
simultaneous  equations  have  given  us  points,  B,  C,  and  Q 
in  Fig.  1.  The  first  two  points  lie  in  the  solution  area  and 
their  coordinates  provide  workable  programs.  Point  Q 
falls  outside  this  solution  area  and  is  unsatisfactory.  But 
that  was  proved  above,  where  we  saw  that  41.1  hours  >  36 
hours. 
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Let  us  next  add  a  third  dimension  to  our  graph.  In  a 
conventional  representation  with  axes  at  90  degrees  from 
each  other,  the  value  of  z  in  a  third  dimension  is  plotted 
for  several  program  points  in  Fig.  2.  Comparing  this  with 
Fig.  1,  we  observe  that  the  z  values  corresponding  to  the 
points  A,  B,  C,  D,  and  0  have  been  shown  as  elevations 


Fig.  2.    Plot  of  z  =  /(Xi,  A2)  in  the  third  dimension  for  Problem  1. 


above  the  plane  established  by  the  coordinate  axes  Xi,  X2. 
These  points,  E,  F,  G,  H,  and  0,  respectively,  have  been 
joined  to  establish  a  bounding  plane.  Or  it  might  be  said 
that  the  plane  EFGHO  is  established  by  the  rays  OE,  OF, 
OG,  and  OH,  of  which  only  the  first  and  last  are  drawn  in 
the  figure.  The  whole  figure  may  be  called  a  convex  poly- 
hedron, and  it  may  be  said  that  the  point  of  maximum 
profit  (i.e.,  the  solution  point)  is  contained  within  this 
polyhedron. 
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The  particular  rays  mentioned  above  are  of  no  im- 
mediate value  to  our  problem,  nor  is  the  plane  OEFGH 
itself  of  primary  interest.  It  is  true  that  it  is  made  up  of 
the  profit  values  corresponding  to  each  workable  program 
and  that  it  shows  point  F  to  represent  the  maximum 


Direction 
■     of 
transla- 
tion 


Fig.  3.     Plot  of  the  z  plane  resting  successively  on  points  H  and  F. 
z=/(Xi,X2). 


profit  whereas  0  represents  the  minimum  profit.  As  such, 
OEFGH  might  be  called  a  "profit  plane."  But  our  par- 
ticular interest  is  the  amount  of  these  profits,  and  the 
maximum  amount  is  given  by  the  perpendicular  BF. 
Again,   if  we  should  become  interested  in  a  less-than- 
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optimal  program,  we  would  wish  to  know  the  length  of 
another  perpendicular.  In  dealing  with  more  compli- 
cated problems  where  more  than  two  X's  require  more 
than  two  dimensions  and,  correspondingly,  a  multi- 
dimensional geometry,  it  is  preferable  to  define  the  term 
"profit  plane"  differently.  Remembering  that  profit  is 
computed  from  our  functional  equation  (/),  we  will  define 
the  term  "profit  plane"  as  a  collection  of  all  points  which 
give  the  same  value  of  z.  And  as  different  programs  are 
successively  tried,  we  will  conceive  of  this  plane  as  being 
translated  parallel  to  itself,  the  length  of  the  perpendicular 
from  the  plane  to  the  origin  being  equal  to  the  profit  at 
any  particular  location  of  the  plane. 

Fig.  3  shows  the  profit  plane  resting  successively  on 
points  H  and  F.  In  this  simple  two-dimensional  example, 
the  perpendicular  to  the  origin  referred  to  above  is,  of 
course,  the  z  axis.  It  may  be  regarded  as  coincidental 
that  all  other  perpendiculars  to  the  base  plane  give  the 
same  value.  In  a  problem  requiring  multidimensional 
representation,  this  condition  would  not  obtain. 

The  Development  of  an  Algorithm 

It  has  been  pointed  out  above  that  the  precise  statement 
of  the  problem  was  in  terms  of  the  inequations 

lOXi  +•   5X2  g  50  (d) 

6X1  +    6X2  g  36  (e) 

4.5Xi  +  18X2  S  81  (g) 

and  the  functional  z  =  f(\u  X2),  or 

9Xx  +    7X2  =  z  (/) 

Now  by  the  simple  device  of  introducing  three  more 
unknowns  which  may  or  may  not  be  equal  to  zero,  we 
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may    convert    these    inequations    to    equations  without 

disturbing  the  precise  statement  of  the  problem.  We  now 
write 

lOXi  +    5X2  +  X3  =  50  (h) 

6X1  +    6X2  +  X4  =  36  (i) 

4.5Xx  +  18X2  +  X5  =  81  (i) 

and  retain  the  same  functional  (/). 

These  equations  are  not  the  equivalent  of  Eqs.  (a),  (b), 
and  (c)  in  that  they  permit  programs  (points  of  solution) 
which  do  not  lie  on  any  of  the  three  lines  of  Fig.  1.  The 
limitation  of  the  program  to  within  the  space  ABC  DO  is 
retained,  but  by  the  adoption  of  these  "slack"  variables 
(X3,  X4,  X5)  the  problem's  requirements  are  more  precisely 
stated.  In  physical  terms,  the  introduction  of  these 
"slack"  unknowns  means  that  we  have  provided  a  method 
of  representing  mathematically  the  failure  to  use  some  of 
the  spare  machine  time  of  any  or  all  of  Lathes  A  and  B 
and  Grinder  G.  Or  we  might  consider  them  to  represent 
imaginary  products  of  zero  profitability.  That  is,  instead 
of  dealing  with  idle  machine  time  per  se,  we  will  make 
products  that  earn  nothing  but  that  also  cost  nothing  to 
make — do  not  wear  out  machines,  use  free  manpower, 
etc. — a  highly  imaginary  type  of  product. 

It  will  be  helpful  in  our  development  to  consider  each 
equation  as  containing  all  of  the  unknowns,  some  with 
coefficients  of  zero,  as  follows: 

10Xx  +    5X2  +  IX3  +  OX4  +  0A5  =  50  (k) 

6XX  +    6X2  +  OX3  +  IX4  +  0X5  =  36  (I) 

4.5Xi  +  18X2  +  OX3  +  OX4  +  IX5  =  81  (m) 

remembering  also 

9Xi  +    7X2  +  0X3  +  OX4  +  0X5  =  z  (fa) 
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This  inclusion  of  zero-value  terms  is  a  device  unknown  to 
common  algebra  and  a  device  to  be  deprecated  there  as 
creating  unnecessary  confusion.  It  is  of  potential  useful- 
ness here  in  making  the  limiting  equations  symmetrical. 

Before  going  further,  since  we  are  developing  a  com- 
putational procedure,  let  us  rearrange  the  order  of  our 
equations  into  a  form  which  will  later  match  that  of  the 
conventional  "simplex."  We  will  place  the  right-hand 
equality  to  the  left,  the  slack  unknowns  next,  and  last, 
the  valued  lambda  at  the  right.  This  arrangement  is 
purely  conventional  and  arbitrary.     Thus: 

'50  =  1X3  +  OX4  +  OX5  +  10Xi  +    5X2      (k) 

Tableau  (a)  \  36  =  0X3  +  1X4  +  0X5  +     6X1  +    6X2       (I) 

81  =  OX3  +  OX4  +  1X5  +  4.5Xi  +  18X2     (m) 

One  program  of  production  which  would  match  this 
Tableau  (a)  or  satisfy  these  spare  machine  time  limitations 
would  be  that  of  no  production  and  would  be  as  shown  in 
Table  1.    We  observe  from  this  table  once  again  that 

TABLE  1 
The  Program  of  Number  of  Lots  That  Matches  Tableau  (a) 


Machine  Hours 

Pin  Type  No. 

1 

2 

Available 

Lathe  A 

0 

0 

50 

Lathe  B 

0 

0 

36 

Grinder  G 

0 

0 

81 

Profit  per  lot 

$9 

$7 

Total  profit 

SO 

$0 

manufacture  of  pin  Type  1  is  more  profitable  than  the 
manufacture  of  any  other  type  of  pin.     In  terms  of  the 
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functional,  this  is  the  same  as  observing  that  the  coefficient 
of  Xi  is  larger  than  that  of  any  other  X  in  the  z  equation. 
Observation  of  this  constitutes  Step  1,  which  is  the  step  of 
deciding  what  type  of  pin  to  "program"  first.  To  put  it 
another  way,  we  decide  to  bring  into  the  program  the 
most  profitable  item. 

We  next  determine  how  many  lots  of  Type  1  pins  can  be 
made.  What,  principally,  is  going  to  limit  the  manu- 
facture of  pin  Type  1?  Practically,  this  is  done  by  dividing 
the  time  available  on  each  machine  by  the  production 
time  required  for  the  manufacture  of  pin  Type  1.  In 
terms  of  our  equations,  this  means  dividing  the  left-hand 
member  by  the  corresponding  coefficient  of  X^  Calling 
this  quotient  6,  we  get 

from  (k)  0i  =  —   =  5 

10 

from  (I)  62  =  —   =  6 

6 

Ol 

from  (m)  03  =  —  =  18 

4.5 

of  which  61  is  the  minimum.  And  this  says  once  again 
that  the  largest  number  of  Type  1  pins  that  can  be  made  is 
five  lots.*  As  Step  2,  then,  a  program  is  chosen  which 
makes  five  lots  of  pin  Type  1.  And  we  use  this  program  to 
modify  the  tableau  by  a  procedure  which  can  at  this  point 
best  be  given  a  pragmatic  justification. 

*  In  particular  instances  in  other  problems  where  negative  coefficients  of 
X  exist,  or  where  negative  "spare  time"  values  are  found  on  the  left  side 
of  the  equations,  6  may  have  a  negative  sign.  Consideration  of  such  a  value 
as  a  limiting  condition  is  useless.  It  would  represent  "negative  production" 
or  the  supplying  of  lots  of  pins  to  the  machines,  running  them  backward, 
securing  raw  material,  and  most  important  of  all — extracting  spare  time ! 
This  may  be  more  clear  if  the  reader  will  consider  the  dimensions  of  6  in  a 
manner  similar  to  the  discussion  following  Tableau  (6)  below. 
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We  ask  next  what  change  would  be  made  in  our  limiting 
conditions  if  a  program  of  making  five  lots  of  Type  1  pins 
were  adopted.  It  appears  that  Tableau  (a),  which  repre- 
sents these  limiting  conditions,  would  have  to  be  modified. 
No  spare  machine  time  would  remain  on  Lathe  A  but, 
rather,  five  lots  of  production  would  fully  use  the  spare 
time  shown  in  Eq.  (k).  The  spare  time  available  on  Lathe 
B  would  be  reduced  to  6  machine  hours.  Similarly  the 
spare  time  on  Grinder  G  will,  by  this  program,  be  reduced 
to  58.5  hours.  The  problem  becomes  one  of  representing 
this  new  set  of  limiting  conditions  in  a  compact  form. 

The  solution  arrangement  which  has  been  adopted  in 
the  simplex  algorithm  does  this  with  a  new  set  of  equa- 
tions identical  in  form  to  the  original  tableau.  The  equa- 
tions replacing  Eqs.  (I)  and  (m)  are  similar  restrictive 
conditions,  except  that  they  are  modified  to  represent  the 
reduced  spare  machine  time  available.  But  the  equation 
replacing  Eq.  (k)  is  quite  different  in  meaning,  if  not  in 
form.  It  represents  lots  of  Type  1  pins  which  have  been 
programmed  and  states  the  restrictions  on  further  program 
possibilities  in  terms  of  these  (five)  lots  rather  than  in  terms 
of  spare  machine  time  available.  Tableau  (6)  defines  the 
constraining  conditions  on  further  programming  in  terms 
of  present  scheduled  lots  of  Type  1  pins  and  remaining 
spare  machine  time.  The  transition  to  this  program  con- 
stitutes Step  2. 

If,  now,  in  the  new  program  we  make  all  we  can  of  pin 
Type  1,  we  will  effectively  divide  Eq.  (k)  by  10  hours  per 
lot  of  Type  1  pins,  as  discussed  above.  We  will  thereby 
develop  a  form  of  equation  which  makes  it  easy  to  solve 
for  X1#  And  we  will  want  a  value  of  Xx  to  substitute  in  the 
other  limiting  equations  in  order  to  move  to  the  situation 
which  describes  the  new  program.    This  gives 

5  =  tV  a3  +  0X4  +  0X5  +  lXj  +  §X2  (*.) 
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or 

Xi  =  5  -  tV As  -  0X4  -  0X5  -  iX2  (kb) 

Once  again  customary  algebraic  manipulation  is  violated, 
and  we  carry  terms  of  zero  value  forward.  The  justifica- 
tion is  that  we  are  making  certain  that  all  operations 
necessary  in  our  final  algorithm  are  recognized,  or  con- 
versely, that  we  have  not  dropped  any  steps  in  a  particular 
solution  because  the  result  of  the  step  was  obviously 
zero. 

Careful  study  of  Eq.  (kb)  may  lead  the  reader  to  wonder 
if  it  is  true.  We  have  said  that  we  propose  to  make  five 
lots  of  Type  1  pins:  thus  Xr  =  5.  But  Eq.  (kb)  shows 
Xi  as  equal  to  five  minus  four  other  quantities.  The  only 
condition,  then,  in  which  Eq.  (kb)  can  be  true  is  the  one  in 
which  the  total  of  these  four  subtracted  quantities  is  zero. 
Recognition  of  this  fact  resolves  the  difficulty.  It  is 
undesirable  to  drop  the  terms  just  because  their  total 
equals  zero.  Certainly  X3  and  X2  might  have  many  values 
which,  when  multiplied  by  — tV  and  —J,  respectively, 
would  give  zero  in  summation.  As  it  turns  out  later,  these 
coefficients  prove  most  helpful  in  determining  succeeding 
programs. 

Substituting  the  value  of  Xx  given  in  Eq.  (kb)  in  Eq.  (I) 
in  order  to  get  a  new  statement  of  condition,  i.e.  a  new 
program,  gives  (in  detail): 

36  =  0X3+1X4+0X5 

+6(5-t1oX3-0X4-0X5-JX2)+6X2     (la) 

36-(6X5)  =  (0-T6o)X3  +  (l-0)X4  +  (0-0)X5 

+0X!  +  (6-3)X2  ft) 

6=  -|X3  +  1X4+0X5+0X1+3X2  ft) 
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Similarly  in  (m) 

81  -0X3  +OX4  +  IX5 

+4.5(5-tVa3-0X4-0X5-|X2)+18X2 

81-  (4.5X5)  =  (0-29o)A3  +  (0-0)X4  +  (l-0)X5 


(-¥) 


+0X!+    18— -   X2  (mb) 


58.5=  -29oX3+0X4  +  lX5+0X1+^X2  (mc) 

Making  these  changes  and  substitutions  constitutes 
Step  2.  Assembling  the  revised  equations  into  a  second 
tableau  similar  to  the  first  yields: 

9    f      5  =  1VX3+0X4+0X5+1X1+JX2  (ka) 

Tableau  (6)  I       6  =  -fXs+lX.+OXs+OXi+SX,  (lc) 

[58.5=  -29oX3+0X4+lX5+0X1+^-3-X2      (mc) 

Let  us  pause  a  moment  and  reconsider  what  we  have 
done  as  far  as  dimensions  are  concerned.  The  value  50 
in  Eq.  (k)  has  the  units  of  hours  time  available.  The 
divisor  used  in  advancing  to  Eq.  (ka)  had  units  of  hours 
per  lot  of  Type  1  pins.  When  the  division  was  made, 
therefore,  the  units  of  the  answer  became  lots  of  Type  1 
pins — the  answer  was,  specifically,  five  lots  of  Type  1  pins. 
If  we  make  five  lots  of  this  type,  we  know  we  will  use 
machine  time  on  Lathe  B  and  Grinder  G.  Looking  first 
at  Lathe  B,  our  original  problem  analysis  indicated  that 
five  lots  of  Type  1  would  use  30  hours,  leaving  6  hours  for 
other  production.  Eq.  (lc)  confirms  this.  Six  hours  of 
machine  time  appear  available  on  Lathe  B.  Similarly 
Grinder  G  will  have  58.5  hours  time  unused,  as  is  seen  from 
Eq.  (mc).  (The  reader  might  well  verify  this  from  the 
original  problem  statement.)     In  the  left-hand  column 
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of  the  group  of  equations  labeled  Tableau  (b),  there 
now  exist  one  figure  (5)  which  has  dimensions  in  lots  of 
Type  1  pins,  and  two  others  which  remain  in  terms  of 
hours.  To  denote  this  distinction,  we  have  placed  the 
symbol  (9)  in  line  with  Eq.  (ka)  on  the  outside  of  the 
bracket.  This  symbol  was  chosen  to  represent  the  profit 
value  of  each  lot  of  this  type  of  pin  brought  into  the  pro- 
gram. If,  now,  we  multiply  a  lot  value  to  the  left  of  the 
equals  sign  but  inside  the  bracket  by  the  lot  profit  value 
in  line  on  the  outside  of  the  bracket  (and  add,  if  there  be 
several  kinds  of  production),  we  will  get  the  profit,  z,  for 
the  particular  tableau.    Thus: 

$9  X  5  =  $45 

This  profit  we  know  to  be  correct.  By  the  adoption 
of  a  trick  of  placement,  we  have  devised  an  automatic  scheme 
for  functional  evaluation,  or  profit  determination,  and  have 
taken  a  small  step  in  the  development  of  our  algorithm. 

Studying  next  the  right  side  of  Eq.  (k),  we  find  that  the 
original  units  were  also  hours — from  the  product  of  lots 
(X)  times  hours-per-lot.  When  the  division  by  10  hours 
per  lot  of  Type  1  pins  was  made,  the  units  of  the  right-hand 
side  of  Eq.  (ka)  also  became  lots  of  Type  1  pins,  making 
the  units  of  (ka)  consistent.  But  the  right  side  of  Eqs.  (lc) 
and  (mc)  remained  in  hours  as  heretofore. 

Tableau  (6)  represents  another  group  of  solutions  to 
the  problem,  but  in  this  tableau  the  constraints  on  any 
solution  have  been  changed.  There  is  no  spare  time  avail- 
able on  Lathe  A.  In  place  of  an  equation  giving  spare 
Lathe- A  time,  an  equation  (ka)  exists  which  calls  for  the 
production  of  five  lots  of  Type  1  pins.  The  other  two 
equations  show  the  reduced  spare  machine  time  that  will 
be  available  if  five  lots  of  Type  1  pins  are  made.  Just  as 
in  the  case  of  Tableau  (a)  we  knew  that  a  program  of 
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making  nothing  would  represent  a  satisfactory  solution 
to  the  constraining  conditions,  we  know  in  this  case  that 
a  program  making  five  lots  of  pin  Type  1  and  no  lots  of 
any  other  type  will  be  satisfactory.  This  is  shown  in 
Table  2.    Since  this  table  shows  some  machine  time  still 

TABLE  2 
The  Program  of  Number  of  Lots  That  Matches  Tableau  (6) 

Machine  Hours 
Pin  Type  No.  1  2  Available 


Lathe  A 

5 

0 

0 

Lathe  B 

5 

0 

6.0 

Grinder  G 

5 

0 

58.5 

Profit  per  lot 

$  9 

$7 

Total  profit 

$45 

$0 

available,  it  is  now  proper  to  inquire  whether  making 
some  lots  of  another  type  of  pin  would  yield  a  greater 
profit.  To  do  this  would  mean  a  reduction  in  the  lots 
of  Type  1  made,  since  our  first  solution  has  used  the  total 
production  time  of  Lathe  A.  Our  problem  is  simple 
enough  so  that  only  a  few  steps  by  cut-and-try  would  give 
the  answer.  Since  there  is  only  one  type  of  pin  left  to 
try,  namely,  Type  2,  we  would  try  a  program  that  made 
one  less  lot  of  Type  1  and  one  lot  of  Type  2  and  check 
for  profit.  Then  we  would  try  a  program  of  two  less  lots 
of  Type  1  and  two  lots  of  Type  2,  and  so  on,  until  from  our 
profit  figures  we  saw  that  we  had  the  problem  bracketed. 
At  this  point  we  might  try  some  fractional  lots  for  greater 
profitability.  However,  we  would  quickly  give  that 
process  up  as  being  too  exhausting  if  the  types  of  pins 
available  for  programming  were  of  the  order  of  four  or 
more.     In  such  a  problem  the  choice  would  have  to  be 
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made  as  to  which  pin  to  experiment  with  first,  and  when 
that  one  had  yielded  its  full  measure  of  profit  increase, 
which  one  to  try  next  and  so  on.  Clearly,  the  problem 
permutes. 

Such  a  choice  being  unnecessary  here,  we  still  try  the  pin 
type  next  most  profitable,  which  is  to  say,  we  repeat  the 
procedure  of  Step  1.  This  step  becomes  a  part  of  our 
algorithm.  When  in  doubt  as  to  which  type  to  program 
next,  choose  the  one  which,  if  it  is  not  yet  in  the  program, 
will  yield  the  highest  profit  in  terms  of  the  functional. 

Taking  Step  2,  we  redetermine  our  productive  limita- 
tions, using  Tableau  (b)  by  a  new  computation  of  0  values. 
As  before,  we  divide  the  left-hand  member  by  the  corre- 
sponding coefficient  of  X2. 

Thus  we  get 

5 


from  Eq.  (ka)* 

0i.i  =  i  =  io 

2 

from  Eq.  (lc) 

^2.1    —    3"    =    2 

from  Eq.  (mc) 

03.1    —    58    X  "63" 

3.71 

from  which  we  observe  that  the  smallest  (limiting) 
quotient  is  2,  or  that  at  the  most,  two  lots  of  Type  2  pins 
can  be  made.  Another  way  of  saying  this  is  to  remark 
that  pin  Type  2  shall  be  brought  into  the  program  in  place 
of  the  limiting  equation  (lc). 

It  is  at  this  point  that  a  very  interesting  bit  of  mathe- 
matical development  occurs  which  may  be  thought  of  as 

*  Dimensions  may  be  of  interest  here.     In  Eq.  (ka)  the  units  of  5  are 
lots  of  Type  1.    The  dimensions  of  \  may  be  computed  as: 

Hours        >        Hours       _  Lot  Type  1 
Lot  Type  2    '   Lot  Type  1  ~  Lot  Type  2 

Hence  B\.\  has  dimensions  of  "lots  Type  2."  In  Eqs.  (Zc)  and  (mc)  the  di- 
mensions are: 

Hours  *  t    4.    ^ 7>  =  Lots  TyPe  2 

Lots  Type  2 
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the  nub  of  the  technique  we  are  trying  to  develop.  We 
desire  to  make  as  many  Type  2  pins  as  the  machine  time 
left  will  permit  us  to  do.  We  have  6  hours  left.  This  will 
permit  two  lots  of  Type  2  pins.  By  bringing  this  type 
into  the  program  in  place  of  the  limiting  equation  (lc),  we 
automatically  reduce  the  quantity  of  Type  1  pins  possible. 
And  we  do  this  by  using  the  same  computational  method 
as  before. 

Repeating  Step  2  and  progressing  now  from  Tableau 
(b),  we  first  divide  Eq.  (lc)  by  the  coefficient  of  X2.  This 
yields 

2  =  -|X3  +  JX4  +  0A5  +  OXi  +  1X2  (ld) 

or 

X2  =  2  +  iX3  -  iX4  +  0X5  +  OXi  (le) 

Substituting  this  value  in  Eq.  (ka)  in  order  to  get  yet 
another  program  gives  (in  detail) : 

5  =  i3o^3+0X4+0X5  +  lX1 

+i(2+ix8-iX4 +0X5+0x0 +ox2    (h) 

5-(|x2)  =(1Lr+1LJ)X3+(0-  J)X4+0X5 

+(1-0)X1+0X2  (kc) 

4=iX3-iX4+0X5  +  lX1+0X2  (kd) 

Also  substituting  in  Eq.  (mc) 

58.5  =  -^X3+0X4+1X5+0X1 

+^(2  +iX3  -  |X4  +0X5  +0X0  +0X2 

{md) 

58.5-(^X2)=[(-^Xi)-^]X3+(0"|X^)X4 

+(1-0)X5+0X1+0X2 

(me) 

27-^Xs— ^X4+1X5+0X! +OX2  (mf) 
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The  revised  tableau  now  is: 


9 
Tableau  (c)     7^ 


4=iX3-iX4+0X5+lX1+0X2  (kd) 

2  =  -  i\3  +iX4  +0X5  +0\!  +  1X2  (ld) 

[27=f£X3-i^X4+lX5+0X1+0X2       (mf) 

The  program  represented  by  this  tableau  of  equations 
calls  for  four  lots  of  pin  Type  1  at  $9  and  two  lots  of  pin 
Type  2  at  $7.  The  profit  from  this  program  will  be  as 
computed  above: 

$9X4  =  $36 
7X2=    14 

$50 

Hence  this  program  is  to  be  preferred. 

Let  us  review  again  the  procedure  we  used  in  deter- 
mining which  pin  type  to  program  first.  We  chose  the 
one  having  the  greatest  profit  value,  namely,  Type  1. 
The  first  choice  always  being  the  type  most  profitable 
from  the  standpoint  of  optimizing  the  functional,  it  is 
always  an  easy  choice  to  make.  In  the  case  of  the  present 
problem,  the  second  choice  was  also  obvious.  There  was 
only  one  other  pin  type  available  for  programming.  But 
if  there  had  been  more  types  available,  the  decision  on 
which  to  program  next  would  have  required  more  study. 
Some  rule  here  would  be  of  help  in  developing  a  system  of 
solution. 

Such  a  rule  can  be  evolved  from  the  process  of  matching 
the  z  equation  resulting  from  each  of  the  tableaux  with  the 
original  z  equation.  To  understand  this  development,  it  is 
better  to  proceed  from  the  second  tableau  rather  than 
from  Tableau  (a),  in  which  no  pins  were  programmed. 
Recall  that  in  Tableau  (b)  we  obtained  the  value  of  $45 
profit  for  z  by  multiplying  the  value  of  five  (lots)  on  the 
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left-hand  side  by  the  profit  value  of  $9.  If  we  now  carry 
the  multiplication  through  on  the  right-hand  side  of  the 
equation,  the  following  expression  for  z  results: 

45  =  t9o  As  +  0X4  +  0A5  +  9Xi  +  |X2  (*) 

This  we  propose  to  compare  with  the  original  equation 
for  z,  item  by  item,  to  see  whether  further  profit  is  possible 
and  which  pin  type  should  be  programmed  to  secure  the 
greatest  return. 

The  comparison  is  made  by  subtracting  each  item  of 
Eq.  (z0)  *  from  each  comparable  item  of  (zb) .    Thus 

45   =  T9o  As   +  0A4  +  0A5   +  9Xi   +  |X2     (zb) 
(-)       0   =     0A3   +OA4   +0A5   +9A!   +7X2     (z0) 
45,  +  T9o A3,  +  0A4,  +  0A5,  +  OAi,  -  |X2 

These  item  values  may  now  be  interpreted  by  straight- 
forward reasoning.  No  change  in  the  programming  of 
pin  Types  1,  4,  or  5  will  increase  or  decrease  the  profit 
of  $45.00.  To  say  this  again:  because  they  have  coeffi- 
cients of  zero,  changing  the  program  to  manufacture  any 
of  Types  1,  4,  or  5  will  not  increase  or  decrease  the  z  value. 
But  with  Types  2  and  3,  a  real  difference  coefficient  exists. 
Making  more  of  Type  3  will  decrease  the  profit  of  $45.00 
by  90  cents  for  each  lot  made.  Conversely,  making  more 
of  Type  2  will  increase  the  profit  of  $45.00  by  $2.50  for 
each  lot  which  is  made.  Hence  a  rule  of  procedure  is 
developed.  Choose  for  the  next  program  to  be  computed 
the  type  of  pin  whose  difference  coefficient  has  the  greatest 
negative  value,  since  this  type  will,  when  programmed,  yield 
the  greatest  increase  in  the  resulting  profit  equation. 

It  may  help  in  seeing  this  if  one  observes  that  a  negative 
value    of    difference  coefficient    means    that    the    profit 

*  Eq.  00)  is  obviously  the  value  of  Eq.  (/a)  when  no  production  is  pro- 
grammed. 
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coefficient  in  the  original  z0  (no-production)  equation  is 
larger  than  the  profit  coefficient  (s)  of  the  programmed 
pin  types. 

Let  us  check  to  see  if  this  rule  would  have  worked  for 
Tableau  (a)  in  advancing  to  (b).  The  subtraction  to  get 
difference  coefficients  would  have  been: 

0   =  0X3   +  0A4   +  OX5   +  OX!   +  0X2      CO 
(-)     0   =  OX3   +  OX4   +  0X5   +  9Xi   +  7X2      GO 
0,       OX3,  +  0X4,  +  0X5,  -  9Xi,  -  7X2 

and  following  our  rule,  we  would  have  chosen  pin  Type  1  as 
the  one  having  the  greatest  negative  value — or  conversely, 
the  greatest  profit  possibility. 

The  rule  has  another  value  which  deserves  mention. 
Since  each  tableau  is  an  iteration  of  a  known  program 
which  satisfies  the  conditions  of  the  problem,  and  since 
it  is  usually  not  our  desire  to  find  all  possible  programs  but 
rather  to  stop  when  a  program  which  optimizes  the  func- 
tional has  been  reached,  let  us  consider  our  rule  in  relation 
to  these  facts.  It  at  once  becomes  obvious  that  if  all 
difference-coefficients  have  become  zero  or  positive  in 
sign,  no  change  in  programming  will  increase  the  left-hand 
side  of  the  equation.  Hence  the  difference  coefficients 
can  also  tell  us  when  an  optimum  program  has  been 
reached.  And  we  have  the  further  rule  that  an  optimum 
program  exists  when  all  difference  coefficients  are  zero  or 
positive. 

Checking  this  rule  against  Tableau  (c),  the  last-devel- 
oped program,  gives: 

36  =  |X3  -|X4  +0X5  -J-9X!  +0X2 
(+)  14  =  -|X3  +fX4  +0X5  +QX!  +7X2 
( = )  50,  +|X3,  +  |X4,  +  0X5,  +  9Xi,  +  7X2 
(-)       0,       +OX3,  +  0X4,  +  0X5,  +  9Xi,  +  7X2 

50,       +IX3,  +  |X4,  +  0X5,  +  OXx,  +  0X2 
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All  coefficients  being  zero  or  positive,  no  type  of  pin  can 
be  introduced  that  will  further  increase  the  profit  of  $50. 
Scheduling  some  lots  of  imaginary  Types  3  and  4  of  pins — 
which  really  means  wasting  more  machine  time — will 
reduce  the  profit. 

Having  thus  found  a  way  of  knowing  that  a  maximum 
program  has  been  reached,  we  turn  our  attention  to  sim- 
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Fig.  4.    Simplex  plot  of  Tableau  (a)  with  computation  of  0  and  za  —  z0. 


plifying  the  computations.  Up  to  this  point  we  have  been 
overcareful  in  preserving  each  detail  of  calculation  in  order 
that  our  final  rules  might  be  checked  against  the  algebraic 
procedure  used  in  development.  It  will  be  convenient  now 
to  recast  each  tableau  in  the  tabular  form  which  has  come 
to  be  known  as  the  "simplex."  Fig.  4  shows  the  original 
all-program  Tableau  (a)  so  arranged. 

A  little  study  shows  this  arrangement  to  be  not  so  very 
different  from  the  forms  of  Eqs.  (k),  (I),  and  (m).  Near 
the  left  is  seen  a  column  denoted  P0  which  lists  the  ma- 
chine hours  available.  The  double  vertical  line  alongside 
may  be  taken  to  denote  the  equals  sign.  The  succeeding 
columns  labeled  P3  .  .  .  P2  correspond  to  the  unknown  X's 
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in  our  equations.  Only  the  coefficients  of  the  unknowns 
are  plotted,  and  where  the  coefficient  is  zero,  the  box 
has  been  left  blank.  At  the  top  of  the  figure,  the  original 
profit  equation  (z0)  has  been  plotted  in  such  a  manner  that 
the  coefficients  fall  above  the  appropriate  X  columns.  The 
original  limiting  equations  have  been  labeled  with  a 
P  designation  corresponding  to  their  distinctive  "slack" 
variables.* 

This  figure  also  shows  a  procedure  for  taking  Step  1  of 
our  method.  Following  the  rule  of  subtracting  the  original 
^-equation  values  from  the  z  equation  determined  by  the 
tableau,  we  observe  the  last  line  to  show  two  profit  possi- 
bilities. Of  the  two,  the  one  having  the  greatest  negative 
value  is  chosen,  and  column  Px  is  investigated  for  limiting 
conditions.  These  computations  of  6  show  at  the  left  of 
the  figure  in  line  with  the  appropriate  equation.  The 
minimum  6  which  describes  the  limiting  productive  condi- 
tion is  also  labeled.  This,  it  will  be  remembered,  was  the 
determining  factor  which  originated  the  equation  manipu- 
lations that  resulted  in  Tableau  (b). 

The  next  step  in  simplification  is  to  state  Tableaux 
(b)  and  (c)  in  a  similar  form.  This  is  done  in  Fig.  5,  which 
also  repeats  Tableau  (a)  and  thus  gives  the  complete 
solution  to  the  problem  in  a  very  compact  form. 

Only  a  few  more  observations  are  needed  to  develop 
the  last  of  the  rules  of  computation  which  will  complete 
the  algorithm  of  the  simplex  tableau.  The  difference 
coefficient  yielding  the  most  negative  value  determines 
which  P  column  shall  be  investigated  for  limiting  condi- 
tions. The  minimum  6  row  is  the  one  chosen  for  all- 
factor  division  by  the  limiting  factor. f    The  value  in  each 

*  Mathematicians  will  recognize  the  identity  matrix. 
t  In  Fig.  5  this  factor  has  been  marked  with  an  asterisk  in  each  of 
Tableaux  (a)  and  (6). 
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box  of  the  minimum  0  row  is  then  divided  by  the  limiting- 
value  (starred)  factor  and  the  results  posted  in  corre- 
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Fig.  5.    Simplex  tableau  solution  of  Problem  1, 


sponding  boxes  of  the  new  tableau.  At  the  same  time, 
any  z  profit  value  appearing  above  the  column  carrying 
the  star  is  posted  to  the  left  of  the  minimum  6  row  in  the 
new  tableau  for  later  use  in  multiplying. 
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With  the  new  row  in  place,  we  proceed  to  compute  new 
values  for  the  rest  of  the  boxes.  The  algorithm  procedure 
is  as  follows: 

1.  Select  a  box  to  be  revalued.    Call  the  old  value  w. 

2.  Move  sideways  in  the  old  tableau  to  the  P  column  to  be 
programmed  in  the  next  tableau  and  note  the  value  in  it. 
Call  it  x. 

3.  Go  to  the  box  in  the  new  tableau  corresponding  to  the  w  box 
of  Rule  1  above.  Then  go  up  or  down  to  the  minimum  6  row 
just  posted.    Note  this  value.    Call  it  y. 

4.  Multiply  x  by  y. 

5.  Subtract  xy  from  w  and  post  this  value  in  the  new  tableau. 

With  some  practice  this  procedure  can  become  very 
rapid.  In  many  cases  it  will  result  in  merely  reposting 
w,  there  being  only  zero  values  in  the  x  or  y  positions. 
But  in  other  cases  it  results  in  revaluation. 

These  particular  rules  were  formulated  after  observing 
the  algebraic  procedure  used  in  the  first  development  of 
Problem  1.  Referring  back  to  Eq.  (lb)  let  us  consider 
the  first  term  after  the  equals  sign.  We  observe  that  the 
corresponding  box  to  be  revalued  is  marked  with  an  (r)  in 
Fig.  5.  The  old  value  in  this  box  was  zero — or  w  =  0. 
Following  Rule  2,  we  move  sideways  to  find  that  x  =  6. 
Observing  Rule  3,  we  find  y  to  equal  tV-  Eq.  (lb)  illus- 
trates Rules  4  and  5,  namely,  the  quantity  within  the 
bracket  is  (0  —  yxr).  And  the  posted  value  corresponds 
to  the  coefficient  — f  shown  in  the  corresponding  position 
of  Eq.  (lc)  or  in  Tableau  (b). 

At  this  point  the  justification  for  carrying  out  the 
algebraic  transactions  of  the  earlier  calculations  in  so 
much  detail  becomes  apparent.  The  reader  is  urged  to 
study  Eqs.  (lb),  (mb),  (kc),  and  (rac),  coefficient  by  coeffi- 
cient, in  the  light  of  the  rules  above,  until  he  is  convinced 
of  the  equivalence  of  the  operations. 
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Problem  2.  Let  the  shop  described  in  Problem  1 
have  the  further  possibility  of  manufacturing  Types 
6  and  7  of  taper  pins,  in  addition  to  Types  1  and  2.  These 
pins  have  the  following  machine  time  and  profit  charac- 
teristics: 


Machine  Time  in  Hours  per  Lot 

Profit  in 

Pin  Type 

Lathe  A     Lathe  B       Grinder  G 

Dollars 

6 

7 

2                2                   1.5 
1                 2                   6.0 

6 
4 

Arrange  this  problem  in  simplex  form  and  solve  for  a 
maximum  profit  program.     (Answer:  $108) 

Problem  3.  Introduce  into  Problem  2  a  management 
policy  requirement  that  at  least  one  lot  each  of  taper  pin 
Types  1  and  2  must  be  scheduled  in  the  spare  machine 
time  available.  Program  the  situation,  using  a  simplex 
tableau  solution  for  maximum  profit.  (Hint:  Does  not 
this  add  the  requirement  Xx  ^  1  and  X2  ^  1?  If  so,  would 
the  equations 

Xi  +  X8  =  1 
and 

X2  +  X9  =  1 

be  added  to  the  three  used  heretofore?) 

Penalty  Coefficients 

If  the  reader  has  taken  the  trouble  to  solve  Problem  3, 
using  the  hint  supplied  above,  he  found  that  the  program 
developed  called  for  18  lots  of  pin  Type  6  and  no  lots  of 
any  other  type.  The  profit  on  this  program  was  identical 
to  that  of  Problem  2.  Moreover  the  program  did  not 
satisfy  the  requirements  of  Problem  3  as  stated,  inasmuch 
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as  no  lots  of  pin  Types  1  and  2  appear  in  the  P0  column. 
Something  is  definitely  wrong  with  the  solution. 

This  problem  was  introduced  to  show  the  necessity  of 
care  in  setting  up  the  mathematical  expression  of  the 
problem  conditions.  There  was  nothing  essentially  wrong 
with  the  statements 

and 


Xi  +  X8  =  1 


x2  +  X9  =  1 

They  merely  were  not  restrictive  enough,  and  the  solution 
bypassed  them  in  maximizing  the  functional.  A  moment's 
reflection  makes  this  obvious.  If  Xi  is  allowed  to  equal 
zero,  then  a  program  of  one  lot  of  X8  will  satisfy  the  first 
equation.  And  similarly,  one  lot  of  X9  will  satisfy  the 
second  equation  with  X2  =  0.  If  the  reader  will  examine 
the  program  shown  in  his  final  tableau  (and  if  he  has 
worked  correctly),  he  will  find  that  it  calls  for  such  "slack" 
production  together  with  the  (18)  very  profitable  lots  of 
pin  Type  6.  Since  this  is  a  more  profitable  program  than 
one  having  some  Type  1  and  Type  2  lots  scheduled,  and 
since  we  did  not  get  the  restrictions  completely  drawn, 
our  algorithm  led  us  to  the  most  profitable  solution,  which 
in  this  case  does  not  fulfill  the  stated  requirements. 

Further  study  shows  no  way  in  which  we  can  readily 
effect  this  degree  of  restriction  using  the  equation  forms 
so  far  developed.  But  the  result  desired  can  be  achieved 
in  another  way.  Suppose  that  some  penalty  were  attached 
to  the  "slack"  variables  X8  and  X9.  It  will  be  remembered 
that  these  "slack"  unknowns  were  introduced  as  imaginary 
products  of  zero  profitability.  If,  now,  scheduling  them 
meant  a  loss  instead  of  a  profit,  the  iterations  of  the 
simplex  tableau  should  readily  eliminate  them  as  various 
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programs  are  developed.    And  this  would  be  particularly 
true  if  the  associated  loss  were  large. 

It  is  not  necessary  even  to  specify  the  amount  of  loss 
to  be  associated  with  an  undesired  variable  which  must 
not  appear  in  a  final  program.  For  our  purposes  let  us 
use  the  letter  M  (for  "mighty")  to  symbolize  a  value 
larger  than  any  other  profit  coefficient,  and  let  us  also 
make  it  negative  in  sign  to  indicate  that  it  is  a  loss  (i.e., 
negative  profit).  Then  the  functional  to  be  associated 
with  Problem  3  would  be 

z  =  9Xj  +  7X2  +  6X6  +  4X7  -  M\8  -  M\9 

(any  other  A's  having  profit  coefficients  of  zero). 

Nothing  new  needs  to  be  added  to  the  simplex  tableau. 
These  negative  profit  values  are  posted  at  the  heads  of 
the  appropriate  columns  and  to  the  left  of  the  appropriate 
rows,  as  was  done  with  other  z  coefficients  heretofore. 
Computation  of  z  for  each  iteration  is  made  in  terms  of 
M's  and  numerical  values.  The  column  to  enter  the  next 
array  as  a  row  is  determined  from  the  greatest  negative 
value  in  the  zt  —  z0  row.  It  is  to  be  expected  that  the 
M-penalized  rows  will  be  replaced  first  and  that  no  such 
row  will  appear  in  the  final  tableau. 

Fig.  6  shows  Problem  3  solved  in  four  iterations  after 
the  introduction  of  this  penalty  device.  It  will  be  observed 
that  some  negative  numerical  values  appear  in  the  last 
zd  —  z0  line.  However,  these  are  subtracted  from  asso- 
ciated, very  large  M  values,  and  hence  the  last  line  is 
entirely  zero  or  positive,  indicating  a  maximal  program. 

Formulation  of  the  Product-Mix  Problem 

As  has  been  previously  observed,  one  of  the  factors 
limiting  the  use  of  linear  programming  as  a  common  en- 
gineering technique   is   the   difficulty  in  expressing  the 
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Fig,  6,    Simplex  tableau  solution  of  Problem  3. 
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problem  situation  in  linear  relationships.  The  problems 
so  far  considered  have  dealt  with  two  types  of  restriction: 
machine  time  available  and  policy  requirements.  Many 
other  types  of  restriction  might  be  added.  One,  which 
has  some  practical  value,  has  a  product-mixture  require- 
ment in  addition  to  the  production  limitation.  The 
development  of  a  mathematical  formulation  of  such  a 
problem  is  illustrated  in  the  following  example.* 

Problem  4.  A  certain  manufacturer  of  screw  fastenings 
found  that  there  was  a  market  for  packages  of  mixed 
screw  sizes.  His  market  research  data  indicated  that  four 
mixtures  of  three  screw  types  (1,  2,  and  3),  properly 
priced,  would  be  most  acceptable  to  the  public.  These 
data  were: 

Selling  Price 
Mixture  Specifications  (^/lb.l 

Not  less  than  40%  Type  1 
A  Not  more  than  20%  Type  2  60 

Any  quantity  Type  3 
Not  less  than  20%  Type  1 

B  Not  more  than  40%  Type  2  25 

Any  quantity  Type  3 
Not  less  than  50%  Type  1 

C  Not  more  than  10%  Type  2  35 

Any  quantity  Type  3 

D  No  restrictions  20 


For  these  screws  the  plant  capacity  and  manufacturing 
costs  are  given  below. 

*  This  type  of  problem  received  early  recognition  in  the  blending  of 
gasolines  in  the  petroleum  industry.  A  simplified  "nut-mix"  example  is 
used  by  Charnes,  Cooper,  and  Henderson  in  their  Introduction  to  Linear 
Programming  (Reference  2)  as  a  medium  for  illustrating  the  simplex 
algorithm. 
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Plant  Capacity  Manufacturing  Cost 

Screw  Type  (lb./day  X  100)  (ff/lb.) 

1  10  50 

2  10  30 

3  6  18 

What  production  shall  this  manufacturer  schedule  for 
greatest  profit,  assuming  that  he  can  sell  all  that  he 
manufactures? 

We  first  state  the  mixture  specification  restrictions. 
Let  A1  indicate  the  pounds  of  Type  1  screws  in  Mixture  A  ; 
similarly,  Bx  the  pounds  of  Type  1  screws  in  B,  etc.  Then 
for  the  mixtures  as  stated  the  restrictions  are: 


i  *u 

B1  feiB 

C,  &|C 

2    &U 

B2  gfB 

^2    =  TO"^ 

As  defined,  the  total  mixtures  will  be: 
A  =  Ai  +.  A2  +  A% 
B  —  Bi  +  B2  +  B% 

C  =  Ci  +  C2  ~\-  Cs 

since  any  quantity  of  Type  3  screws  may  be  used  to  com- 
plete the  mixtures. 

We  next  substitute  the  last  three  definitional  relation- 
ships in  the  restrictive  inequations,  giving: 

A1  g  Ui  +  U»  +  Ms  or  0  &  -Ux  +  U»  +  Uz 

M  S  Ui  +  U2  +  U3  or  0  g  tAt  -  iA2  +  iA3 

Bt  S  m  +  iB2  +  |B3  or  0  S  -iBt  +  \B2  +  |B3 

B2  S  f B,  +  |B2  +  f  B3  or  0  g  IB,  -  f B2  +  f  B3 

Cx  ^  JC,  +  §C2  +  iC3  or  0  ^  -JCi  +  |C2  +  1C3 

C  ^  x ow  +  To ^2  +  To ^3  or  0^  to ^i  ~~  To ^2  +  T0C3 
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By  multiplying  through  where  appropriate  by  a  —  1,  the 
inequalities  may  all  be  turned  in  the  same  direction: 

o  a  -Ui  +U2+U3 
0  a  -iAi  +U2-U3 
0  g  -#B,  +  iB2  +  iB3 
0  &  -*Bi+lB2-lB3 

0    ^        TOW   T"  To^2         T¥^3 

These  six  inequations  state  the  restrictions  in  so  far  as 
mixtures  are  concerned.  We  next  write  the  production 
limitations.    These  are: 

Screw  Type  1  =  Ax  +  Bx  +  d  +  Dx  S  10 

Screw  Type  2  =  A2 +B2  +  C2  +  D2  £  10 

Screw  Type  3  =  A3  +  B3  +  C3  +  D3  ^    6 

We  now  rewrite  our  inequations  in  the  customary  X 
terminology  and  introduce  " slack"  unknowns  to  reduce 
the  inequations  to  equations.*    This  gives: 

-fAi  +  f  X2  +  IX3  +  x13  =  0 

—  6^X1    +   5^X2    —    5"X3    +  X14    =    0 

-*x4+lx5+ix6  +  x15  =  0 

—  5^4  +  5"X5  —  f  Xe  +  X16  =  0 

-ix7  +  |x8  +  *x9+x17  =  0 

—  To  ^7  +  T0^8  —  To"X9  +  Xis   =  0 

*  The  order  of  X  substitution  for  letter  notation  used  here  was  alphabeti- 
cal by  subscripts,  with  the  higher  subscript  values  being  reserved  for  the 
slack  variables.  A  tableau  form  thus  evolves,  which  is  similar  to  the  ones 
so  far  used.    This  order  is  a  matter  of  convenience  only. 


THE  REPLACEMENT  ROW  DIFFICULTY  37 

Ai  +  X4  +  X7  +  Xio  +  x19  =  10 

X2  +  X5  +  X8  +  Xn  +  X20  =  10 

X3  +  X6  +  X9  +  X12  +  X2i  =  6 

Our  equations  are  now  in  a  form  for  ready  arrangement 
into  a  simplex  tableau.  However,  we  still  must  formulate 
a  functional — the  relationship  to  be  optimized.  Observing 
that  selling  price  and  manufacturing  cost  are  given,  we 
write  a  profit  functional: 

z  =  selling  price  of  mixtures  minus  the  cost  of  components 
or 

z  =  60  (Xx  +  X2  +  X3)  +  25(X4  +  X5  +  X6) 

+  35  (X7  +  X8  +  X9)  +  20(X10  +  Xn  +  X12) 
-50(Xi  +X4+X7+X10) 
-30(X2  +X5+X8+Xn) 
-  18  (X3  +X6+X9+X12) 

This  equation  simplifies  to 

z  =  10XX  +  30X2  +  42X3  -  25X4  -  5X5 
+  7X6  -  15X7  +  5X8 
+  17X9  -  30X10  -  10Xn  +  2X12 

Fig.  7  shows  these  equations  arranged  in  the  regular 
simplex  form  with  the  trivial  no-production,  no-profit 
program.  Fig.  8  displays  an  optimum  program  in  terms 
of  the  functional. 

The  Replacement  Row  Difficulty 

When  one  is  working  with  the  simplex  algorithm,  a 
situation  occasionally  arises  in  which  no  row  satisfies  the 
requirement  that  6  be  positive  and  either  zero  or  greater. 
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Such  is  the  case  in  the  tableau  presented  in  Fig.  9.  Here 
the  only  column  showing  a  negative  value  in  the  z  —  z0 
row  is  the  one  headed  P6.  In  this  column  the  coefficients 
are  all  zero  or  negative.  Dividing  these  coefficients  into 
the  Pq  " limiting  conditions"  column  yields  values  of  6 
which  are,  as  shown,  either  infinite,  indeterminate,  or 
negative.  Clearly  an  impasse  exists.  No  row  will  satisfy 
the  requirements  of  our  method. 
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Fig.  9.    Illustration  of  the  replacement  row  difficulty. 


This  situation  results  from  one  of  two  conditions.  The 
first  (and  in  the  author's  experience,  most  common)  reason 
for  this  situation's  developing  is  that  a  mistake  has  been 
made  in  computations.  Unfortunately,  such  mistakes  are 
usually  several  tableaux  back  and  are  often  hard  to  find. 
Usually  it  is  best  to  revert  to  the  basic  tableau  and  begin 
afresh,  stopping  at  the  end  of  each  section  of  the  work  for 
comparison  with  the  first  computation.  Although  it 
appears  longer,  this  technique  has  often  proved  to  be  the 
most  economical. 

If  there  has  been  no  error  in  computation,  the  reason 
for  this  impasse  must  lie  in  the  problem  itself.     It  is 
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possible  that  the  restrictive  conditions  as  set  forth  in  the 
original  linear  equations  did  not  define  a  meaningful  pro- 
gramming problem.  This  most  often  occurs  when  the 
requirements  have  been  so  stated  that  optimization  of  the 
functional  would  approach  infinity.  In  terms  of  a  profit 
functional  such  as  is  used  in  economic  problems,  this  means 
that  the  restrictive  equations  permit  unlimited  activity 
in  one  or  more  profitable  directions  and  that  there  are 
not,  perhaps,  enough  limiting  conditions  stated  in  the 
original  problem. 

The  simplex  tableau  cannot,  of  itself,  point  out  the  lack 
in  the  initial  problem  statement.  Nothing  thus  far  de- 
veloped can  take  the  place  of  the  mind  of  the  investigator 
who  states  the  situation.  But  by  failing  to  be  workable 
at  some  particular  point,  the  technique  does  bring  to  his 
attention  the  fact  that  he  has  stated  a  situation  which 
cannot  be  optimized.  Further  study  (perhaps  with  the 
aid  of  the  tableaux)  may  enable  him  to  see  why  the 
problem  is  empty. 

The  general  rule  may  therefore  be  stated:  If  all  the 
coefficients  in  the  replacing  column  are  zero  or  negative  and 
no  computational  error  has  been  made,  the  original  statement 
of  conditions  did  not  define  a  meaningful  linear  program- 
ming problem. 

Degeneracy 

In  advancing  to  a  new  tableau  of  equal  or  greater 
profitability,  the  simplex  algorithm  requires  that  the 
values  in  the  P0  column  be  inspected  for  the  next  imme- 
diately limiting  condition.  This  is  done  by  computing  a 
6  value  for  each  limiting  equation  (row)  and  selecting  the 
lowest  positive  value  (negative  production  not  permitted 
— see  footnote,  page  15).    But  on  many  occasions,  "ties" 
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for  minimum  6  values  occur  between  two  or  more  equa- 
tions. Such  a  situation  is  called  "degeneracy"  and  indi- 
cates that  our  solution  method  has  failed  to  be  automatic 
in  its  application. 

When  such  a  situation  exists,  two  or  more  of  the  rows 
appear  equally  well  qualified  to  be  replaced.  If,  in  their 
P0  column,  they  have  zero  values  (as  is  the  case  when 
the  6  tie  is  zero),  it  can  be  assumed  that  they  do  not  have 
immediate  programming  relevance.  The  trouble  with 
ignoring  them  on  this  basis  (as  was  often  done  in  the  past) 
is  that  better  programs  may  not  be  independent  of  such 
equations  and  that  the  road  to  better  programs  may 
involve  their  use,  albeit  at  no  immediate  functional  in- 
crease. The  problem,  then,  is  to  find  a  rule  of  procedure 
to  guide  us  in  our  choice  of  the  row  to  replace  and  thus 
make  the  method  automatic. 

Dr.  Charnes*  has  provided  such  a  rule,  and  in  so  doing 
has  eliminated  one  of  the  major  objections  to  the  simplex 
method.  Considering  the  value  of  P0  in  each  row-equa- 
tion to  be  made  up  of  the  sum  of  the  products  (coefficient 
multiplied  by  X)  from  left  to  right  across  the  tableau,  he 
suggests  that  in  the  case  of  6  ties  for  Po's,  one  compare 
further  the  individual  component  products  and  resolve 
the  tie  on  the  basis  of  these  products,  selecting  the  smallest 
while  proceeding  from  left  to  right.  An  example  should 
help  to  make  this  clear. 

In  Fig.  10  the  Tableau  (e)  requires  that  the  replacing 
column  be  P8.  Theta  values  as  shown  at  the  left  reveal  a 
zero  tie  between  rows  PX7  and  P18.  Following  Charnes's 
rule,  the  degeneracy  has  been  handled  by  the  comparison 
shown  at  the  bottom  of  the  tableau.  Here,  in  addition  to 
dividing  the  total  P0  element  by  the  incoming  coefficient 

*  Reference  2. 
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which  resulted  in  equivalence  of  limitation,  each  compo- 
nent product  making  up  the  P0  has  been  similarly  checked 
for  equivalence.  Proceeding  from  left  to  right,  the  tie 
was  first  broken  in  column  P17  and  at  this  point,  since 
row  Pi8  is  the  smaller,  it  was  chosen  as  the  row  to  be 
replaced.  The  principal  reason  for  moving  from  left  to 
right  is  to  insure  that  the  "slack"  coefficients  shall  be 
inspected  for  non-equivalence  first. 

The  rule  is  general.  A  tie  cannot  obtain  in  every 
column  across  the  tableau  because  each  equation,  by 
virtue  of  defining  a  specific  and  unique  limiting  condition, 
cannot  be  identical  with  another  equation.  This  is  defined 
in  the  original  "basic"  square  of  the  original  "slack" 
coefficient  setup.  At  the  intersection  of  its  own  column, 
each  original  equation  has  the  coefficient  of  one.  All  other 
elements  in  that  column  must  be  zero.  And  this  is  where 
the  tie  is  usually  broken. 

Programming  with  Exact  Restrictions 

All  of  the  problems  examined  thus  far  have  had  their 
original  requirements  expressed  in  the  form  of  inequa- 
tions. If  the  simplex  technique  were  limited  to  such 
situations  only,  its  value  as  an  analysis  technique  would 
be  narrow  indeed.  Fortunately,  it  is  possible,  by  a  simple 
device,  to  extend  it  to  cover  situations  in  which  some  exact 
requirements,  defined  by  equations,  must  be  satisfied. 
Let  us  consider  such  a  situation.  We  choose  one  similar 
to  Problem  2. 

Problem  5.  A  certain  shop  having  machine  time 
available  wished  to  schedule,  most  profitably,  four  types 
of  tapered  fastening  pins.  Machine  tool  time  available 
and  profit  values  are  given  below. 
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Lot  Manufacturing 

Time  (Hours)  for 

Pin  Type 


Machine  Type 


Lathe 
A 


Lathe 
B 


Grinder 
G 


Profit 

in 
Dollars 
per  Lot 


1 
2 

6 

7 

Machine  time  (hours) 
available 


10 
5 
2 

1 


50 


36 


4.5 

18.0 

1.5 

6.0 


81.0 


It  is  further  required  that  all  of  the  time  available  on 
the  Grinder  G  be  used,  since  this  machine  cannot  be  used 
for  production  of  any  other  kind,  whereas  the  lathes  may 
be  used  in  making  screws  (perhaps).  What  program  will 
satisfy  these  requirements  most  profitably? 

Stating  the  requirements  mathematically,   as  hereto- 
fore, gives 

10Xi  +    5X2  +     2X6  +    X7  ^  50 
6X1  +    6X2  +     2X6  +  2X7  ^  36 


4.5Xi  +  18X2  +  1.5X6  +  6X7  =  81 

the  functional  to  be  maximized  being 

9Xi  +  7X2  +  2X6  +  4X7  =  z0 

Our  conditions  provide  two  inequations  and  one  equation. 
It  is,  of  course,  the  latter  which  we  have  not  yet  learned 
to  handle. 

To  make  equations  amenable  to  this  technique,  recourse 
is  had  to  the  device  of  introducing  "artificial  slack"  vari- 
ables with  a  penalty  coefficient.  The  penalty  coefficient 
insures  that  such  "artificial  slack"  is  rapidly  removed  as 
the  iterations  of  the  solution  proceed.    We  have  already 


46 


LINEAR  PROGRAMMING 


used  the  penalty  coefficient  device  for  removing  unwanted 
rows  in  the  final  tableau  of  Problem  3;  hence  no  new 
operational  procedure  is  introduced.  For  this  problem, 
the  tableau  equations  become 

50  =  X3  +   lOXj  +    5X2  +     2X6  +    X7 

36  =  X4  +     6X1  +    6X2  +     2X6  +  2X7 

81  =  X5  +  4.5Xi  +  18X2  +  1.5X6  +  6X7 

The  functional  to  be  used  is 

z  =  -M\5  +  9Xi  +  7X2  +  2X6  +  4X7 

The  first  array  for  this  problem  is  displayed  in  Fig.  11. 
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Fig.  11.     Tableau  (a)  for  Problem  5. 


Carrying  the  iterations  through  as  before  yields  a  pro- 
gram as  shown  in  Table  3.  If  the  reader  will  carry  the 
simplex  method  through,  he  will  see  that  the  penalized  row 
is  the  first  to  be  replaced  and  that  the  tableaux  develop, 
with  each  succeeding  one  showing  an  increase  in  the  profit 
functional.  With  this  penalty  coefficient  device  available, 
problems  having  some  exact  equivalence  restrictions  and 
those  with  broader  limitations  are  equally  amenable  to  the 
simplex  technique. 
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TABLE  3 

An  Optimum  Program  for  Problem  5 


Machine  Time  Used 

Pin             Lots 
Type       Scheduled 

Lathe 
A 

Lathe 
B 

Grinder 
G 

Profit 
($) 

1  2 

2  0 

6  0 

7  12 

Total  (hours) 
Hours  unused 

20 
0 
0 

12 

32 

18 

12 
0 
0 

24 

36 
0 

9 
0 
0 

72 

81 
0 

18 
0 
0 

48 

($)66 

Development  of  Alternate  Programs 

There  are  often  instances  arising  in  which  an  optimum 
program  is  developed  (as  shown  by  the  Zi  —  z0  row  being 
all  zero  or  positive)  without  all  of  the  columns  in  the 
array  having  been  used.  One  properly  wonders,  then, 
if  other  optimum  programs  might  not  exist  were  other 
columns  substituted  for  appropriate  rows.  It  might 
well  be  that  a  different  optimum  program  could  be  found, 
which  would  be  preferable  to  the  one  first  obtained  for 
reasons  not  stated  in  the  original  problem  conditions. 

Careful  analysis  can  help  one  determine  if  such  alternate 
optimum  programs  are  possible.  We  observe  first  that 
the  alternate  program (s)  must  have  a  zt  value*  equal  to 
the  z  value  of  the  optimum  program.  Briefly,  if  an 
alternate  program  is  also  optimum,  zovt.  =  Zi,  or  Zi  — 
zopt.  =  0-  This  observation  suggests  that  we  consider 
more  closely  the  relationship  between  the  successive 
z  values  as  tableaux  iterate  in  the  direction  of  maximiza- 

*  The  subscript  i  indicates  that  the  z  value  is  computed  from  any  tab- 
leau i. 
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tion.    Since  this  can  be  done  most  easily  by  example,  let 
us  consider  Tableaux  (b)  and  (c)  of  Fig.  5  (page  28). 

In  Tableau  (c),  zc  =  50;  in  (b),  zb  =  45.  How  was  the 
increase  computed?  Observing  that  0  equaled  2,  we  see 
that  the  following  steps  were  taken.  The  value  of  P0  in 
row  Pi  was  recomputed  according  to  rule 

5  -  (J  X  9)  =  4 

In  the  second  row,  P0  equaled  0,  and  reflection  tells  us  that 
this  is  always  true  for  the  row  which  is  being  replaced  by 
a  column.  In  computing  the  value  of  zc,  4  was  multiplied 
by  the  value  9  in  column  z0,  and  0  was  multiplied  by  a 
corresponding  value  of  7.  These  together  totaled  50  =  zc. 
We  now  indulge  in  a  little  algebraic  manipulation.  The 
increase  in  z  as  we  progressed  from  zb  to  zc  may  be  written 

zc-zb  =  [(9  X  4)  +  (7  X  2)]  -  (9  X  5) 

Substituting  4  =  5  —  ^0,  and  also  2  =  0,  from  above, 

zc-Z6  =  9X5- fl9  +  7<9 -9X5 

zc-  zb  =  7d  -  10 

zc  -  zb  =  0  (7  -  |) 

which  gives  a  statement  of  the  increase  in  terms  of  0 
and  two  numerical  values. 

If  we  could  generalize  these  numerical  values,  a  general 
statement  of  the  z  increase  effected  by  replacing  a  row 
with  a  column  would  result.  This  is  not  too  hard  to  do. 
The  value  7  is  the  z0  profit  coefficient  of  the  column  that 
is  being  introduced  as  a  row.  Let  us  designate  such  a 
column  as  Pk*  The  statement  of  the  increase  might  then 
be  written  (using  a  consistent  subscript  notation) 

zc  —  zb  =  6(zok  —  1) 

*  We're  being  cute,  k  is  for  "kum  in."  The  o  is  to  indicate  that  it's  the 
original  profit  value  of  the  "kum  in"  column.    Coincidentally,  ok. 
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Observe,  now,  the  zk  value  of  the  Pk  column  in  the 
illustration  (Fig.  5)  from  which  we  have  been  generalizing. 
It  agrees  in  value  with  the  remaining  numerical  value  in 
the  equation  just  above.  And,  indeed,  it  should,  for  the 
f  in  the  equation  above  came  from  the  value  of  \  in  the 
Pk  column  multiplied  by  the  corresponding  profit  coeffi- 
cient in  the  z0  column.  This  was  not  a  singular  coinci- 
dence. Had  there  been  more  than  just  one  replacement 
row  (Pi)  in  Tableau  (b),  all  others  would  have  had  their 
values  in  the  P0  column  modified  (as  the  5  was  modified) 
by  6  and  by  the  corresponding  value  (like  J)  in  the  Pk 
column.  And  these  modified  values  would  have  been 
multiplied  by  their  corresponding  profit  coefficient  in  the 
z0  column  before  addition  to  secure  a  value  (like  f )  in 
the  Zk  box. 

Hence  we  generalize  the  statement  of  the  z  increase  as 

Z(i+D  —  Zi  =  6(zok  —  zk) 

where  zok  is  the  original  profit  value  of  a  replacing  column, 
zk  is  the  sum  of  the  profits  in  the  Pk  column  of  the  tableau, 
and  6  is  a  ratio  of  the  P0  value  of  a  row  being  replaced  to 
the  replacing  column. 

And  so  a  useful  relationship  has  been  discovered.  The 
value  in  the  Pk  box  of  the  zt-  row  when  subtracted  from  the 
original  profit  value  of  the  Pk  column  will,  after  being 
multiplied  by  6,  determine  the  profit  increase  of  the  tableau 
iteration.  But  we  observe  further  that  this  is  one  of  the 
values  which  is  always  computed  as  the  last  line  of  a 
tableau  (the  Zi  —  z0  line)  in  order  to  determine  whether 
an  optimum  program  has  been  reached,  or  to  determine 
in  what  direction  to  advance.  Only,  in  these  routine 
calculations,  the  direction  of  subtraction  is  reversed,  and 
so  the  sign  of  the  answer  is  reversed. 
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Let  us  now  return  to  our  original  purpose  in  analyzing 
for  a  statement  of  tableau  increase.  We  had  observed 
that  once  an  optimum  tableau  was  reached,  the  possi- 
bility of  alternate  optimum  programs  depended  upon 
the  increase 

%>i  ~  2opt.    —   0 

But  any  increase  between  tableaux  is  required  by  our 
algorithm  to  be  equal  to 

6(zok  —  zk) 

Hence  the  condition  for  alternate  optimum  programs  is 
that 

%i  —  Sopt.  =  0  =  6(zok  —  Zk) 

And  this  condition  can  obtain  either  when  6  is  zero  or  when 
the  statement  within  the  brackets  is  zero. 

Thus  we  may  readily  inspect  any  optimum  tableau  for 
the  possibility  of  alternate  equivalents.  We  inspect  first 
the  columns  which  have  not  been  entered  as  rows.  If  any 
of  these  have  in  their  z{  —  z0  box  a  zero  value,  the  condi- 
tion above  is  satisfied  and  it  should  be  possible  to  use  this 
column  as  a  Pk  row  to  effect  another  optimum  program. 
Again,  we  may  review  the  unused  columns  to  find  if  they 
can  in  any  instance  yield  a  6  value  of  zero.  If  so,  the  zero- 
increase  condition  of  the  equation  above  is  satisfied,  and 
another  optimum  tableau  is  available. 

But  alternate  tableaux  effected  as  a  result  of  zero  6 
values  do  not  achieve  new  usable  programs.  The  only 
way  in  which  a  zero  6  becomes  possible  is  when  the  row 
value  under  PQ  is  zero.  Thus,  although  a  new  column  is 
entered  into  the  row  equations,  it  comes  in  with  a  P0  value 
of  zero  and  does  not  change  any  other  P0  value.  Since  the 
chosen  optimum  program  is  taken  from  the  P0  column, 
the  effectiveness  of  zero  6  values  in  determining  alternate 
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optimum  programs  is  illusory.  Nevertheless  this  does  not 
invalidate  the  mathematics  above  nor  the  effectiveness 
of  the  (zi  —  z0)  row  as  an  alternate  tableau  selection 
device. 
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Fig.  12.     Two  tableaux  representing  alternate  optimum  programs. 


Fig.  12  illustrates  the  use  of  the  (Zi  —  z0)  device  in  the 
selection  of  an  alternate  optimum  program.  In  Tableau 
(/)  the  first  maximum  program  has  evolved.  All  values 
in  the  (zf  —  z0)  row  are  zero  or  positive.  Inspecting  this 
tableau  for  alternate  optimum  program  possibilities  shows, 
first,  that  columns  P4,  P7,  and  P8  have  not  been  entered  as 
rows  or  programmed.  Of  these,  only  column  P8  has  a 
zf  —  z0  value  of  zero.*  Hence  it  is  the  only  possibility 
for  an  alternate.  In  array  (g),  it  has  replaced  row  Pu. 
As  was  anticipated,  the  functional  remained  constant  at 


*  Remember  that  (z/  —  z0)  = 
B{zi  -  z0). 


(zok  —  Zk).    In  general,  -6{zok  —  zk)  = 
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600  units.  Inspection  for  6  values  from  columns  P4  and  P7 
yields  no  6  values  equal  to  zero.  It  may  safely  be  con- 
cluded that  for  the  problem  represented  by  Fig.  12,  only 
two  optimum  programs  are  available. 

Problem  6.  Analyze  Problems  4  and  5  for  alternate 
optimum  programs. 

Problem  7.  A  certain  production  scheduling  problem 
involved  two  machines  which  had  180  and  90  hours,  re- 
spectively, available.  Three  products,  A,  B,  and  C,  were 
to  be  made,  which  involved  machine  times  as  shown  in 
the  top  left  box  of  Fig.  13.  In  this  box  the  X's  represent  the 
unknown  program.  Some  other  restrictions  were  also 
involved. 

This  problem  was  maximized  in  six  arrays.  Tableau  (/) 
is  also  shown  in  Fig.  13.  The  bottom  left  box  shows  the 
corresponding  maximum  machine  allocation  program. 
What  alternate  maximum  programs  are  available? 

Analysis  for  Non-Optimum  Programs 

The  development  of  an  optimum  program  by  the  use 
of  the  simplex  algorithm,  together  with  the  development 
of  the  relationship  explained  in  the  section  on  Alternate 
Programs,  provides  a  rapid  means  of  investigating  the 
relative  value  of  non-optimum  programs.  In  many 
situations  it  may  be  as  desirable  to  know  the  amount  of 
"loss"  that  will  result  from  a  less-than-optimum  program 
as  it  is  to  know  the  value  of  the  functional  given  by  the 
optimum  program.  Once  an  optimum  program  has  been 
evolved  and  the  value  of  the  functional  for  the  optimum 
condition  computed,  a  "peak"  exists  to  which  the  "profit- 
ability" of  less-than-optimum  programs  may  be  compared. 
It  may  well  be  that  if  the  "loss"  resulting  from  a  particular 
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less-than-optimum  program  is  not  too  great,  it  will  be 
preferred,  for  some  reason (s)  not  stated  in  the  original 
array,  to  the  optimum  program.  And  since  a  rapid  means 
of  analysis  for  this  "loss"  figure  exists,  it  is  obvious  that 
the  usefulness  of  the  simplex  algorithm  is  greater  than 
has  heretofore  been  evident. 

It  is  first  observed  that  in  the  evolution  of  the  optimum 
program  a  number  of  inferior  programs  (profitwise)  are 
iterated.  These  may  be  readily  compared  for  the  amount 
of  "loss"  resulting  from  their  use  by  subtracting  their  Zi 
value  from  the  2opt.  value.  Other  programs  not  yet  stated 
must,  of  course,  come  from  the  introduction  of  unused 
columns  into  the  row-equation  position  to  form  new  arrays. 

The  development  of  the  expression  for  an  increase  in 
functional  value  in  the  section  on  Alternate  Programs 
yielded 

Increase  =  6{zok  —  zk)  =  —  0(z*  —  z0). 

The  expression  to  the  right  is  the  one  of  present  interest. 
It  states  the  two  terms  which  govern  the  increase  to  be 
expected  when  a  particular  column  is  used  to  replace  a 
particular  row.  For  the  simplex  array,  the  bracketed 
term  (Zi  —  z0)  is  computed  at  the  bottom  of  the  particular 
column.  (If  it  is  negative,  an  increase  may  result.)  Theta 
is  the  ratio  we  get  by  dividing  the  P0  value  of  the  to-be- 
replaced  row  by  the  intersecting  row-value  of  the  partic- 
ular column  under  study.  Hence  the  increase  in  z  which 
could  be  effected  by  an  iteration  can  be  known  even  before 
the  computation  is  made.  By  using  2opt.  for  ziy  an  opti- 
mum "bench  mark"  value  is  set,  to  which  all  other  pro- 
grams may  be  compared  profitwise  without  the  necessity 
of  making  the  tableau  iterations  (except  as  one  might 
wish  for  a  check  in  a  few  cases). 
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Let  us  use  Fig.  12  to  illustrate  this  device.  In  Tableau 
(g)  an  optimum  program  exists,  as  is  seen  from  the  zg  —  z0 
line  values,  all  being  positive  or  zero.  The  columns 
headed  P4  and  P7  are  not  in  the  array.  Suppose  we  ask, 
"What  loss  would  result  from  introducing  P4  into  a 
program?'7 

Note  that  for  column  P4 

(^opt.    —   %o)    =    \Zg  —  2o)    =    13 

and  that  the  6  values  resulting  from  replacing  various  rows 
by  this  column  would  be  as  shown  in  Table  4. 

TABLE  4 

Increase  in  Functional  Value  Obtained  from  Programming 
Column  P4  in  Tableau  (g)  of  Fig.  12 


Increase 

Row 

e 

(Zff  -   Zo) 

-  (Zg   -    Zo)0 

Ps 

100 

13 

-1,300 

p2 

100 

13 

-1,300 

p6 

0 

13 

0 

p5 

0 

13 

0 

Pi 

100 

13 

-1,300 

Ps 

331 

13 

-433i 

The  products  resulting  from  various  substitutions  are  also 
shown  in  this  table,  and  since  they  all  represent  reductions 
(negative  increases)  from  the  zopt.  value  of  600,  it  is  seen 
that  P8  represents  the  only  replacement  row  leaving  any 
profit.  Subtraction  of  the  increase  of  433 J  from  the 
optimum  shows  that  a  program  with  column  P4  replacing 
row  P8  will  give  a  z  value  of  166f .  Thus  an  analysis  has 
been  made  without  an  iteration. 

In  Fig.  14,  the  tableau  matching  this  program  has  been 
partially  developed  to  the  point  where  it  can  serve  as  a 
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check.    The  remainder  of  the  array  is  left  unfinished  as 
an  exercise  for  the  reader. 
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Fig.  14.     Check  of  a  non-optimum  program.    (Array  (h)  is  left  unfinished 
as  an  exercise  for  the  reader.) 

Problem  8.     Using  the  data  of  Problem  7,  determine 
the  second-  and  third- best  programs.    Answer  =  260. 


Further  Illustration  of  Increase  Formula  Usage 

Another  illustration  of  the  use  of  this  method  of  rapid 
analysis  for  non-optimum  programs  may  be  found  by 
referring  back  to  Problem  3,  page  30.  Solution  of  this 
problem  without  the  use  of  the  penalty  coefficient  device 
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led  to  an  optimum  solution  which  did  not  contain  the 
restrictive  conditions 

Xi  ^  1     and     X2  ^  1 

in  the  final  array.    This  tableau  is  shown  in  Fig.  15. 
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Fig.  15.     An  improperly  restricted  optimum  array  for  Problem  3. 


Let  us  assume  that  we  had  developed  such  a  solution 
and  were  unaware  of  the  penalty  coefficient  device.  By 
the  use  of  the  increase  relationship,  we  could  determine 
relatively  quickly  how  best  to  force  into  the  program 
one  lot  each  of  Type  1  and  Type  2  pins.  Since  zf  —  z0  for 
column  Px  is  smaller  than  for  column  P2,  we  first  analyze 
row  by  row  for  its  replacement,  as  shown  in  Fig.  16. 
Clearly  the  minimum  loss  ($9)  will  result  if  this  column 
is  substituted  for  row  P8.  So  doing  will  set  up  a  new 
array,  as  shown.  It  has  been  bracketed  and  called  (g). 
Because  we  require  the  final  program  to  contain  P2,  we 
carry  the  computations  only  so  far  as  necessary,  compute 
columns  P0  and  P2,  and  leave  the  array  in  skeleton  form. 
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Fig.  16.    Use  of  increase  formula  for  rapid  solution. 
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Zg  plus  increase  equals  Zh  or  $88. 

Analyzing  as  before  for  the  minimum  loss  resulting 
from  introducing  column  P2,  we  see  that  it  occurs  when 
row  P9  is  replaced.     And  adding  this   negative  profit 
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to  our  previously  determined  value  gives  a  total  profit 
of  $88  for  this  program. 

Problem  9.  Two  products,  A  and  B,  each  require 
three  sequential  departmental  production  operations  in 
the  course  of  manufacture.  The  first  two  operations  are 
performed  in  production  departments  W  and  X.  The  last 
may  be  done  in  departments  Y  and  Z  or  by  overtime  in  Z 
(which  is  designated  Zf  below).  Hence  there  are  three 
schedules  of  manufacture  possible  for  each  product,  as 
shown  in  Table  5  below.  This  table  gives  the  unit  pro- 
duction times  for  each  department,  the  unit  profit,  and 
the  time  available  for  production  in  each  department. 
It  is  required  that  the  most  profitable  schedule  of  manu- 
facture be  determined  for  products  A  and  B. 

TABLE  5 

Unit  Process  Times  in  Minutes  per  Piece  for 
Problems  9  to  14,  Inclusive 


Product  A 

Product  B 

Produc- 

Produc- 

tion 

tion 

Sched- 

Sched- 

Sched- 

Sched- 

Sched- 

Sched- 

Hours 

Depart- 

ule 

ule 

ule 

ule 

ule 

ule 

Avail- 

ment 

I 

II 

III 

I 

II 

III 

able 

W 

6 

6 

6 

4 

4 

4 

34,000 

X 

5 

5 

5 

2 

2 

2 

18,000 

Y 

8 

3 

10,000 

Z 

7 

3 

3,500 

Z' 

10 

4 

14,000 

Unit 

Profit 

$1.20 

$1.05 

$0.95 

$0.65 

$0.45 

$0.50 

Problem  10.  Methods-improvement  work  resulted  in 
the  production  time  (Table  5)  for  Product  A  in  Depart- 
ment X  being  reduced  to  4  minutes  per  piece,  and  in 
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Department  Z  to  6  minutes  per  piece.     What  optimum 
schedule  now  exists  and  what  profit  results? 

Problem  11.  If  a  methods  change  in  the  manufacture 
of  Product  B  reduces  production  time  in  Department  Y 
to  2  minutes  per  piece,  what  optimum  schedule  results 
at  what  profit? 

Problem  12.  Assume  that  competition  has  resulted  in 
each  unit  profit  shown  in  Table  5  being  cut  by  $0.07.  What 
optimum  schedule  and  what  profit  would  now  obtain? 

Problem  13.  Let  4,000  hours  more  overtime  be  avail- 
able in  addition  to  that  shown  in  Table  5.  What  is  the 
optimum  schedule  for  these  products  and  at  what  profit? 

Problem  14.  Let  it  be  further  required  that  exactly 
100,000  units  of  Product  A  be  made.  Let  other  conditions 
remain  as  stated  in  Table  5.  What  schedule  should  be 
set  for  maximum  profit? 

Minimization 

All  of  the  illustrations  used  and  all  of  the  problems 
introduced  thus  far  have  dealt  with  the  case  of  maximizing 
a  functional  equation.  In  all  instances  this  functional 
has  been  a  money  (profit)  value.  In  most  of  the  writing, 
however,  we  have  been  careful  to  speak  of  this  process  as  one 
of  optimization  rather  than  maximization,  thus  implying 
that  the  functional  might  also  be  minimized.  With  the  pos- 
sibility of  problems  in  minimization  being  also  amenable 
to  this  technique,  its  value  Avould  be  further  enhanced. 
Such  is  indeed  the  case.    We  consider  this  topic  now. 

In  passing,  we  should  further  observe  that  in  the  devel- 
opment of  our  algorithm  through  the  use  of  an  example 
which  used  profitability  in  money  as  a  criterion,  we  did  not 
thereby  restrict  our  method  only  to  such  profit  conditions. 
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The  restriction  placed  upon  our  functional  was  not  one  of 
units  but  rather  that  it  be  expressable  in  equation  form 
and  use  the  same  unknowns  that  were  used  in  the  restric- 
tive equations.  Our  next  illustrative  example  will  use  a 
functional  whose  units  are  not  money  but  time.  As  before, 
an  example  will  be  used  to  introduce  the  use  of  the  simplex 
tableau  in  a  situation  requiring  the  minimization  of  a 
functional.  And,  as  before,  the  problem  will  be  made  so 
simple  as  to  render  it  readily  soluble  by  "cut-and-try." 

Problem  15.  A  certain  machine  shop  had,  during  a 
given  month,  spare  machine  time  available  on  three  lathes 
as  shown  in  Table  6  below.  Four  articles  could  be  sched- 
uled to  these  machines,  and  these  articles  had  production 
requirements  for  this  given  month,  as  the  table  shows. 
Machining  times  in  hours  have  been  entered  in  the  body 
of  the  table.  Minimum  machine  use  was  desired.  What 
program  of  machine  assignments  should  the  shop  use?* 

TABLE  6 
Data  for  Problem  15 


Product 

ion 

Time  (hr 

./pc.) 

Available 

Article 

I 

II 

III 

IV 

Machine  Hours 

Lathe  A 

3 

3 

2 

5 

80 

Lathe  B 

4 

1 

1 

2 

30 

Lathe  C 

2 

2 

3 

1 

130 

Required 

Production 

10 

40 

50 

20 

*  This  problem  can  be  solved  very  quickly  by  a  straightforward  "di- 
minishing-returns"  method,  when  the  data  are  arranged  in  table  form  and 
the  condition  to  be  minimized  is  so  simple  a  function  of  the  requirement. 
One  would  first  schedule  all  of  Article  IV  on  Lathe  C,  all  of  Article  I  on 
Lathe  C,  etc.,  striving  always  to  use  a  minimum  of  machine  time.  It  will, 
however,  because  it  can  be  so  easily  checked,  still  serve  as  an  introduction 
to  minimization. 
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We  first  set  up  a  table  of  unknown  production  quantities 
in  the  same  form  as  the  table  of  data  above  to  help  in 
keeping  our  unknowns  orderly.  Such  a  table,  given  below, 
becomes,  once  the  unknowns  are  determined,  a  statement 
of  the  optimum  program  of  machine  assignments. 

TABLE  7 
Unknown  Program  for  Problem  15 


Production  Pro 

gram 

(Pieces) 

Available 

Article 

I 

II 

III 

IV 

Machine  Hours 

Lathe  A 

Xi 

x2 

x3 

x4 

80 

Lathe  B 

x5 

x6 

x7 

x8 

30 

Lathe  C 

x9 

Xio 

Xn 

Xi2 

130 

Required 

Production 

10 

40 

50 

20 

Reference  to  the  problem  statement  and  the  tables 
permits  the  following  statement  of  requirements  in  the 
form  of  inequations  and  equations: 

3Xx  +  3X2    +  2X3    +  5X4    ^    80 

4X5  +    X6    +    X7    +  2X8    <    30 

2X9  +  2X10  +  3Xn  +    X12  ^  130 

\i  +   x5  +  x9  =    10 

X2    +    X6    +    X10  =    40 

X3    +    X7    +    Xn  =    50 

X4    +    X8    +    X12  =    20 

Adding  the  slack  unknowns  gives: 

3Xx  +  3X2    +  2X3    +  5X4    +  X13  =    80 

4X5  +    X6    +    X7    +  2X8    +  X14  =    30 
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2X9  +  2Xio  +  3Xn  +   x12  +  x15  = 

130 

Xi      +      X5     +      X9     +  ^16    = 

10 

X2      +      X6      +      X10    +  X17    = 

40 
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X3    +    X7    +    Xn  +  Xis  =    50 
X4   +    X8   +    X12  +  X19  =    20 

The  requirement  stated  in  the  functional  is  to  be  one 
of  minimum  machine  use.  This  says  that  the  sum  of 
Xi  +  X2  +  •  •  •  +  X12  scheduled,  when  multiplied  by  their 
corresponding  machine  time  values,  shall  be  as  small 
as  possible.  Also,  since  we  have  four  exact  restrictions 
in  our  tableau  of  equations,  we  wish  to  insure  that  X16, 
X17,  Xi8,  and  X19  have  values  of  zero  or  do  not  appear  in  the 
final  array.  As  before,  we  achieve  this  by  giving  them  a 
penalty  coefficient  in  the  functional  equation.  But  since 
we  desire  minimization,  the  penalty  coefficient  in  this 
case  will  be  a  large,  positive,  but  unspecified,  machine 
time  value.  As  before,  we  will  use  the  letter  M  to  repre- 
sent this  value.    The  functional  to  be  minimized  thus  is: 

z0  =  3Xx  +  3X2  +  2X3  +  5X4  +  4X5  +  X6  +  X7  +  2X8  + 

2X9  +  2X10  +  3XU  +  Xia  +  M\1Q  +  M\17  +  MX18  +  MX19 

The  first  tableau  statement  of  this  problem  is  shown  in 
array  (a)  of  Fig.  17.  As  it  stands,  it  is  identical  in  form 
to  those  tableaux  we  have  studied  heretofore.  The  first 
change  in  our  algorithm  develops  when  we  proceed  to  the 
second  array.  This  time  we  wish  to  progress  in  a  direction 
which  will  minimize  the  functional,  and  by  the  most  direct 
route.  Essentially  this  means  that  we  will  program  first 
the  column  which  causes  the  least  usage  of  machine  time. 
Inspection  of  the  z0  line  in  Tableau  (a)  shows  that  either 
column  P6,  P7,  or  P12  could  be  chosen,  since  they  tie  at  a 
minimum  value  of  1  hour  per  piece  machine  time.     In 
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terms  of  the  line  za  —  z0  which  we  have  become  used  to 
inspecting,  this  is  equivalent  to  advancing  on  the  basis 
of  the  most  positive  difference  coefficient.  And  this,  being 
the  reverse  of  our  previous  maximization  idea,  gives,  intu- 
itively, a  feeling  of  correctness. 
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Fig.  18.     Two  minimum  programs  for  Problem  15. 

To  advance  to  Tableau  (b)  we  have  arbitrarily  chosen 
column  Pi2.  Inasmuch  as  negative  production  is  still 
inadmissible,  we  will  continue  to  use  the  minimum  positive 
6  as  our  criterion  in  determining  the  row  to  be  replaced. 
Hence  the  selection  is  row  P19  with  a  6  value  of  20.  Pro- 
gressing from  array  to  array,  then,  as  in  the  original 
algorithm,  but  choosing  the  columns  to  replace  rows  upon 
the  basis  of  the  maximum  positive  Zi  —  z0,  will  effect  the 
minimization  of  this  problem  in  six  steps,  giving  either  of 
the  two  equivalent  programs  (/)  or  (g)  shown  in  Fig.  18. 
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(The  intervening  steps  which  led  to  this  minimum  of  190 
machine  hours  are  omitted  for  a  reader  exercise.) 

It  is  known  that  minimization  has  been  reached  when 
the  Zi  —  z0  line  has  all  values  zero  or  negative,  since  this 
permits  no  possibility  of  a  further  reduction  by  the  intro- 
duction of  any  other  column  into  the  program. 

As  before,  an  alternate  optimum  program  exists  when 
any  column  not  in  the  array  has  a  difference  coefficient  of 
zero.  Transition  from  array  (/)  to  (g)  in  Fig.  18  was  made 
on  this  basis.  Similarly,  any  increase  over  the  minimum 
machine  hours  given  by  less-than-optimal  programs  can 
be  determined  from  the  product  of  the  Zi  —  z0  value  under 
that  column  multiplied  by  the  6  value  of  the  replaced  row. 

Problem  16.  Develop  the  arrays  missing  between 
Tableaux  (a)  and  (/)  given  in  Figs.  17  and  18. 

Problem  17.  Using  Tableau  (g),  determine  quickly  the 
smallest  increase  in  machine  time  required  if  column  Pi 
is  entered  in  the  final  array. 

TABLE  8 
Data  for  Problem  18 


Product 

ion  Time  (hr./pc.) 

Available 
Machine 

Production  Order 

I 

II 

Ill 

IV 

V 

VI 

Hours 

Machine  No.  1 
Machine  No.  2 
Machine  No.  3 
Required  pro- 
duction, pieces 

3 
4 
2 

10 

3 

1 
2 

40 

2 
1 
5 

60 

5 
2 
1 

20 

2 

2 
1 

20 

1 

1 
2 

30 

80 

30 

160 

270 

Problem  18.     In  the  body  of  Table  8  above,  the  pro- 
duction time  in  hours  per  piece  is  given  for  each  of  six 
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production  orders  which  require  operations  on  three 
machines.  At  the  right  of  the  table,  the  available  machine 
hours  are  listed.  At  the  bottom  of  the  table  the  size  (in 
pieces)  of  each  production  order  is  given.  Schedule  this 
production  so  as  to  minimize  the  machining  time. 

The  Problem  of  Distribution  Costs 

One  area  to  which  linear  programming  was  early  and 
fruitfully  applied  was  that  of  minimizing  distribution 
costs.  In  form,  these  problems  are  identical  to  the  illus- 
tration (Problem  15)  used  to  introduce  minimization. 
They  differ  only  in  units  used.  Instead  of  minimizing 
hours  of  machine  time,  they  minimize  the  cost  of  trans- 
portation from  one  location  to  another.  Quite  often 
these  problems  have  involved  so  few  items  that  they  have 
been  readily  soluble  (as  was  the  case  with  Problem  15)  by 
cut-and-try.  In  other  cases,  the  variables  have  been 
numerous  enough  to  warrant  an  attack  with  the  simplex 
technique. 

We  have  included  such  a  problem  immediately  below, 
which  may  be  handled  in  10  to  30  minutes  by  cut-and-try 
or  by  a  diminishing-returns  procedure,  albeit  without  any 
guarantee  that  a  true  minimum  has  been  reached.  It  is, 
however,  detailed  enough  to  deserve  analysis  with  the 
simplex  algorithm. 

Problem  19.  A  certain  company  making  a  bulk  product 
operated  five  manufacturing  plants  well  located  from  the 
standpoint  of  raw  materials,  water,  power,  and  other 
production  factors,  but  relatively  remote  from  their  mar- 
kets. In  order  to  serve  its  customers  better,  the  company 
had  established  three  warehouses  economically  located 
with  respect  to  principal  markets.  The  problem  which  the 
company  faced  was  how  much  product  from  each  plant 
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to  send  to  each  warehouse  in  order  that  shipping  costs 
might  be  minimized. 

Study  yielded  the  summarized  data  shown  in  Table  9 
below  with  regard  to  plant  capacities,  market  require- 
ments (both  in  tons  per  day),  and  shipping  costs  (in 
cents  per  ton). 


TA] 

Data  for 

3LE  9 

Problem  19 

Unit  Shipping  Cost  from 

Market 
Require- 
ments 

Plant  to  Warehouse 

Plant 

I 

II 

III 

IV 

V 

Warehouse  A 

25 

33 

51 

16 

40 

60 

Warehouse  B 

40 

25 

33 

19 

20 

160 

Warehouse  C 

60 

36 

26 

39 

17 

55 

Plant 

Capacity 

75 

92 

40 

10 

68 

It  was  obvious  that  the  warehouses  could  not  accept  or 
dispose  of  more  than  the  given  daily  tonnage  (on  the 
average),  but  the  plants  could  be  operated  at  any  desired 
level  below  capacity.  How  many  tons  should  be  sent  from 
each  plant  to  each  warehouse  in  order  to  achieve  the  lowest 
over-all  shipping  cost? 

Reduction  of  this  data  to  linear  equations  is  identical 
in  form  to  that  shown  in  Problem  15,  and  it  will  not  be 
repeated  here.  The  functional  will  have  three  lambda 
values  with  M  coefficients,  in  order  to  insure  the  removal 
of  the  "artificial  slack"  variables  introduced  in  the  case 
of  the  equations  representing  the  exact  restriction  on  ware- 
house receipts.  The  functional  will  be,  of  course,  the  sum 
of  the  products  of  the  unknown  shipping  quantities  multi- 
plied by  their  respective  shipping  costs.     If  we  assign 
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lambda  variables  in  "book  reading"  order  in  the  data 
table,  reserving  the  later  serial  values  for  the  slack  vari- 
ables, the  functional  equation  to  be  minimized  would  be 

z0  =  25XX  +  33X2  +  51X3  +  I6X4  +  40X5 

+  40X6  +  25X7  +  33X8  +  19X9  +  20X10 
+  60Xn  +  36X12  +  26X13  +  39X14  +  17X15 
+  M\16  +  M\17  +  MX18. 

Tableau  (a)  for  this  problem  is  shown  in  Fig.  19.  The 
final  minimum  shipping  cost  program  is  given  in  Table  10. 
Iteration  of  the  intervening  tableaux  is  left  to  the  reader 
as  an  exercise.  Since  all  lambda  coefficients  are  valued 
at  one,  these  arrays  can  be  developed  very  rapidly.  The 
final,  minimal,  array  will  be  recognized  as  one  having 
all  Zi  —  z0  values  as  zero  or  negative  in  sign.  Inspection 
of  this  array  according  to  our  rules  will  show  that  only 
one  minimal  program  exists  for  this  problem. 


TABLE  10 

The  Minimum 

-Cost 

Shippin 

g  Pro 

3RAM 

in  Tons 

per  Day 

from  Plant  to  Warehouse 

FOR 

Problem 

19 

Market  Re- 

Plant 

I 

II 

III 

IV 

V 

quirements 

Warehouse  A 

60 

60 

Warehouse  B 

5 

92 

10 

53 

160 

Warehouse  C 

40 

15 

55 

Plant  Capacity 

75 

92 

40 

10 

68 

Unused 

10 

0 

0 

0 

0 

Further  Studies 

As  was  indicated  in  the  introduction  to  this  book,  the 
simplex  algorithm  as  a  mathematical  device  is  now  approx- 
imately five  years  old.    The  first  public  announcement  of 
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its  application  to  industrial  problems  was  made  about  four 
years  ago.  Because  of  this  and  because  of  the  advanced 
mathematical  derivation,  the  number  of  different  types  of 
industrial  situations  in  which  it  has  been  found  applicable 
is  not  large,  although  the  situations  are  in  themselves 
consequential.  And,  as  has  been  previously  mentioned, 
it  is  anticipated  that  a  general  widespread  understanding 
of  the  method  will  result  in  many  more  novel  applications. 

In  addition  to  the  problems  already  discussed,  which  in 
our  opinion  best  illustrate  the  types  most  often  studied 
in  industry  to  date,  some  others  have  been  reported  which 
illustrate  unique  applications  of  this  or  similar  techniques. 
One  of  these  is  an  algorithm  developed  especially  for,  and 
limited  to,  the  " transportation"  problem  that  was  dis- 
cussed above  under  the  heading  of  Distribution  Costs. 
It  will  be  remembered  that  the  solution  of  the  problem 
of  distribution  costs  by  the  simplex  technique  developed 
very  rapidly  because  all  lambda  coefficients  were  valued 
at  unity.  Because  the  restrictive  equations  had  this 
characteristic,  another  algorithm  requiring  very  little 
computation  has  been  developed  for  problems  of  this  type. 
It  has  been  named  the  " stepping  stone  method."  Once 
the  matrix  array  has  been  set  from  the  limiting  equations, 
this  technique,  like  the  simplex,  can  be  carried  out  by 
clerical  help  following  a  set  of  rules.  An  explanation  of 
the  "  stepping  stone  method"  will  be  found  in  Reference  6. 

One  of  the  most  interesting  applications  of  the  simplex 
technique  (and  one  which  constitutes  a  contribution  to 
problem  formulation)  to  what  has  been  shown  above  was 
recently  made  by  Charnes,  Cooper,  and  Ferguson*  in  their 
study  of  executive  compensation.  For  this  study,  data 
were  assumed  to  exist  in  the  form  of  ratings  of  company 
executives  on  a  variety  of  listed  management  traits.    The 

*  Reference  5. 
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compensation  of  these  executives  was  also  known.  It 
was  desired  that  a  scheme  be  evolved  whereby  one  might 
know  what  portion  of  an  executive's  pay  was  being  granted 
for  each  trait ;  that  is,  whether  there  existed  in  the  data, 
as  a  result  of  past  company  practice,  some  implicit  and 
precise  but  unrecognized  relationship  between  total  pay 
and  the  relative  amount  of  certain  desirable  traits  pos- 
sessed by  individual  executives. 

The  authors  formulate  this  problem  first  in  nonlinear 
mathematical  terms  and  then  reduce  the  problem,  by  the 
use  of  additional  equations,  to  a  linear  form  which  is 
amenable  to  a  simplex  solution.  From  the  optimum 
program  thus  developed  they  are  readily  able  to  write 
an  estimating  equation  which  serves  to  determine  the 
value  the  company  places  upon  the  degree  of  possession 
of  a  given  management  trait.  The  estimating  equation 
could  be  used  to  set  the  compensation  for  a  new  (and 
similar)  executive  job,  provided  only  that  the  job  be 
describable  in  terms  of  the  amount  of  each  of  the  given 
list  of  management  traits  required  on  the  part  of  its 
holder.  As  such,  the  study  represents  a  real  contribution 
to  the  art  of  job  evaluation  and  wage  administration. 

The  article  referred  to  also  explains  the  mathematical 
theory  involved  in  the  reduction  of  the  original  nonlinear 
problem  to  linear  form.  The  existence  of  this  reduction 
possibility  suggests  that  ingenuity  can  extend  the  simplex 
technique  to  a  still  wider  variety  of  management  problems. 
Charnes  et  al.  feel  that  such  extension  is  especially  likely 
in  the  area  of  market  analysis — perhaps  with  relation  to 
the  setting  of  sales  quotas. 

It  will  be  recalled  that  one  of  the  first  management 
areas  to  which  the  simplex  algorithm  was  applied  was 
located  in  the  refining  industry.     Symonds'  recent  book* 

*  Reference  8. 


FURTHER  STUDIES  73 

presents  a  collection  of  the  types  of  problems  in  this 
industry  to  which  the  method  is  applicable.  They  include: 
problems  of  product  mixes  considering  fuel  oils,  lubrica- 
tion oils,  and  aviation  gasolines  when  given  particular 
crude  oil  capacities;  problems  of  production  rates  and 
inventory  vs.  seasonal  requirements ;  transportation  prob- 
lems, including  crude  oil  allocation,  so  that  optimizing 
involves  both  maximizing  and  minimizing;  and  some 
applications  of  game  theory  as  formulations  to  be  solved 
by  linear  programming.  The  student  who  wishes  to 
learn  more  about  the  successful  application  of  this 
algorithm  to  industrial  problems  will  find  this  book  of 
considerable  interest.  With  the  background  he  has  ob- 
tained from  the  volume  in  hand,  he  should  have  little 
difficulty  in  following  Symonds7  developments. 

A  survey  concerning  the  present  status  of  linear  pro- 
gramming and  one  which  relates  it  to  the  future  will  be 
found  in  Dantzig's  paper  titled  "Recent  Advances  in 
Linear  Programming."*  In  this  article,  the  inventor  of 
the  simplex  tableau  considers  the  application  of  linear 
programming  in  several  general  areas.  One  reported  is 
where  capabilities  do  not  greatly  exceed  requirements, 
as  was  the  case  in  most  examples  used  in  the  present 
text.  A  second  section  of  the  article  deals  with  the 
application  of  linear  programming  to  other  combinatorial 
types  of  problems.  Herein  he  feels  that  the  methods 
developed  for  the  " assignment"  problem  and  the  " travel- 
ing salesman"!  problem  are  reasonably  adequate  and 
that  further  progress  will  be  in  the  direction  of  applying 

*  Reference  7. 

f  The  "traveling  salesman"  problem  is  to  find  the  best  route  for  visiting  n 
cities  so  as  to  minimize  the  total  distance  covered.  The  "assignment" 
problem  is  a  special  case  of  the  transportation  problem  and  is  one  of  op- 
timally assigning  m  men  to  n  jobs,  knowing  the  value  of  having  the  ^th  man 
on  the^th  job. 


74  LINEAR  PROGRAMMING 

linear  programming  methods  to  new  areas.  It  is  suggested 
that  there  may  well  be  important  applications  still  to  be 
made  in  chemistry,  economics,  engineering,  and  other 
fields  of  mathematics — applications  of  a  type  not  yet 
invented. 

Dantzig  considers  the  last  situation  perhaps  the  one 
of  greatest  interest.  The  size  of  the  problem  matrix 
grows  rapidly  as  restraining  conditions  are  added.  This 
may  make  solution  impossible  timewise  even  when  the 
most  modern  electronic  equipment  is  available.  It  is 
reported  that  a  general-purpose  simplex  code  (for  elec- 
tronic computers)  now  exists  which  can  handle  a  matrix 
of  the  order  of  100  equations  and  1,500  unknowns  with 
machine  times  varying  from  5  to  8  hours.  However,  a 
matrix  of  this  size  may  still  be  too  small  (or  the  time  too 
long)  for  some  business  and  management  problems.  Hence 
Dantzig  logically  suggests  that  future  developments  must 
be  in  the  direction  of  simplifying  the  mathematical  models, 
taking  advantage  of  any  special  properties  a  problem  may 
have,  and  using  any  mathematical  shortcuts  that  ingenuity 
can  devise  before  applying  the  simplex  procedure.  This 
suggestion  implies  that  the  more  complicated  problems 
are  often  going  to  require  the  services  of  a  mathematician 
to  get  them  reduced  to  manageable  size. 

It  is  hoped  that  the  reader,  as  a  result  of  his  experience 
with  this  text,  will  be  encouraged  to  study  more  formal 
mathematical  presentations  of  this  subject — particularly 
since  he  has  now  achieved  some  comprehension  of  a  new 
method.  However,  if  this  presentation  has  left  him  with 
the  feeling  that  further  work  should  be  restricted  to  those 
who  have  a  greater  mathematical  aptitude,  we  still  trust 
that  he  has  achieved  some  new  skill  in  recognizing  the 
types  of  management  and  business  problems  that  are 
amenable  to  mathematical  formulation.     And  we  hope 


REFERENCES  75 

that  he  is  now  better  prepared  to  present  them  to  the 
mathematician  for  expert  assistance. 
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